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FOREWORD 


When  the  National  Bureau  of  Standards  was  established  more  than  80  years  ago, 
it  was  given  the  specific  mission  of  aiding  manufacturing  and  commerce.    Today,  NBS 
remains  the  only  Federal  laboratory  with  this  explicit  goal  of  serving  U.S.  industry 
and  science.    Our  mission  takes  on  special  significance  now  as  the  country  is  respond- 
ing to  serious  challenges  to  its  industry  and  manufacturing—challenges  which  call 
for  government  to  pool  its  scientific  and  technical  resources  with  industry  and 
universities. 

The  links  between  NBS  staff  members  and  our  industrial  colleagues  have  always 
been  strong.    Publication  of  this  new  Industrial  Measurement  Series,  aimed  at  those 
responsible  for  measurement  in  industry,  represents  a  strengthening  of  these  ties. 

The  concept  for  the  series  stems  from  the  joint  efforts  of  the  National  Confer- 
ence of  Standards  Laboratories  and  NBS.    Each  volume  will  be  prepared  jointly  by  an 
industrial  specialist  and  a  member  of  the  NBS  staff.    Each  volume  will  be  written 
within  a  framework  of  industrial  relevance  and  need. 

This  publication,  A  Primer  for  Mass  Metrology,  represents  the  first  of  what  we 
anticipate  will  be  a  long  series  of  collaborative  ventures  that  will  aid  both 
industry  and  NBS. 


Ernest  Ambler,  Director 
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1.  INTRODUCTION 

Much  of  the  mass  program  at  the  National  Bureau  of  Standards  has  been  documented 
during  the  last  two  decades  [1-8].    What  we  attempt  here  is  to  explain  the  relevance 
of  these  publications  to  mass  metrology  conducted  at  laboratories  whose  standards  are 
calibrated  by  NBS.    We  try  to  begin  with  basic  physics  concepts  deriving  the  essential 
results  with  a  minimum  of  rigor.    The  reader  should  always  consult  the  appropriate 
reference  for  a  more  sophisticated  treatment. 

We  have  also  tried  to  use  a  consistent  notation.  Consequently,  our  notation  may 
vary  slightly  from  that  found  in  the  references  cited  above. 

Many  examples  are  provided  in  the  belief  that  understanding  of  a  general  case  is 
much  easier  if  a  special  case  is  understood  well.  Cautionary  statements  are  provided 
at  places  in  the  text  where  important  practical  complications  may  exist. 

In  general,  we  have  strived  to  be  understandable  rather  than  scholarly  where  the 
combination  was  elusive. 

CAUTION:    This  document  is  incomplete.    Many  extremely  important  topics—manu- 
facture, specification,  cleaning,  handling,  and  storage  of  weights  being  among 
them—have  been  omitted.    The  first  two  of  these  subjects  are  treated  in  ref.  [14]. 

Our  ultimate  goal  is  to  outline  the  way  in  which  mass  measurements  can  be  made 
to  acceptable  metrological  accuracies.    There  are  three  major  concepts  which  must  be 
explored  but  which  are  somewhat  interrelated  -  what  one  means  by  "mass;"  how  one  can 
use  a  measuring  device  to  determine  the  mass  of  an  unknown  object  in  terms  of  an 
accepted  unit,  for  example,  the  kilogram;  and  how  to  assign  a  realistic  uncertainty 
to  the  results.    It  will  be  useful  to  introduce  some  basic  equations  of  physics  in 
order  to  make  our  arguments  precise.    The  first  major  step  is  to  explore  the  basic 
properties  of  a  mass-measuring  device,  or  balance,  and  then  see  how  the  device  can  be 
put  to  the  service  of  metrology.    Mathematically,  this  means  that  the  first  step  is 
to  derive  the  basic  equations  of  weighing.    The  remainder  of  our  efforts  modify  the 
basic  equations  in  order  to  make  them  practical.    We  also  show  their  relevance  to 
mass  metrology. 
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2.    BASIC  MASS  EQUATIONS 


In  order  to  derive  the  basic  equations  of  weighing,  we  need  a  precise  notion  of 
mass.  For  the  purpose  of  this  discussion  we  will  introduce  mass  through  the  familiar 
equation: 

f  =  Ml 

(Force  =  Mass  x  Acceleration)  (1) 

That  is,  each  object  possesses  a  property  called  "mass"  which  appears  in  eq.  (1)  as 
the  constant  of  proportionality  between  a  force  applied  to  that  object  and  the  result- 
ing acceleration  of  the  object.    Note  that  both  force  and  acceleration  are  vector 
quantities:    they  have  a  direction  associated  with  them.    One  should  also  note  that 
mass  is  always  a  positive  number  (that  is,  for  example,  the  acceleration  is  always  in 
the  same  direction  as  the  gravitational  force,  never  opposite  to  it).    These  comments 
about  mass  are  consistent  with  the  qualitative  idea  of  the  mass  of  an  object  being  a 
measure  of  the  "amount  of  substance"  in  the  object. 

Unfortunately,  our  intuitive  notions  of  "mass"  are  often  confused  with  "weight." 
Such  confusion  is  unacceptable  for  science  and  metrology.    Therefore  we  will  now  see 
how  the  notion  of  weight  differs  both  logically  and  practically  from  that  of  mass. 

We  can  take  it  as  an  experimental  fact  that  over  a  small  plane  area  of  the 
earth's  surface  (|uch  as  the  space  occupied  by  a  metrology  laboratory)  the  accelera- 
tion of  gravity,  g,  is  essentially  constant.  Since  g  is  a  vector  quantity,  constant 
§  implies  that,  over  a  small  area  of  the  earth's  surface,  the  direction  of  §  is  also 
practically  constant.  That  is,  the  acceleration  vectors  are  all  parallel  and  define 
the  direction  "down."  In  this  approximation,  we  can  replace  §  by  a  numerical  con- 
stant, g. 

The  gravitational  force  on  an  object  of  mass  M  is  then: 

Fj  =  Mg    .  (2) 

The  weight  of  an  object  of  mass  M  is  defined  as  F, .  That  is,  weight  is  a  force,  not 
a  mass. 

Using  W  =  F^ ,  we  have 

W  =  Mg    .  (3) 

From  what  was  said  about  g,  it  is  clear  that  weight  is  not  a  constant  property 
of  matter,  but  depends  on  location.    Consider,  for  example,  a  body  of  mass  M  taken 
out  into  deep  space  such  that  the  gravitational  forces  are  negligible  (they  are  never 
zero).    In  this  case,  we  can  approximate  (3)  as 

W  =  Mg  0 

which  holds  since 

g     0;  (M  f  0)  . 


This  definition  was  adopted  for  international  use  by  the  General  Conference  for 
Weights  and  Measures  (C.G.P.M.)  in  1901. 
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Hence,  the  body  still  has  mass  M  but  its  weight  is  zero.    It  is  only  because  most 
people  live  exclusively  on  the  surface  of  the  earth  that  weight  and  mass  are  often 
confused. 

Let  us  now  consider  the  action  of  a  scale  or  balance  -  the  device  which  we  will 
use  to  measure  mass.    There  are  a  great  many  schemes  in  use  for  constructing  scales 
and  balances.    All  have  their  strengths  and  weaknesses  but  most  are  designed  to 
approximate  the  workings  of  the  following  hypothetical  device: 

This  device  is  a  "black  box"  which  has  a  pan  on  which  to  place  an  object  and  an 
indicator  which  reads  zero  when  the  pan  is  empty  and  some  other  number  when  an 
object  is  placed  on  the  pan  (fig.  1).    Our  ideal  scale  has  the  following  two  de- 
sirable properties: 

1)  The  indicator  always  reads  the  same  amount  to  any  desired  precision  when 
the  same  object  is  put  on  the  pan  at  the  same  conditions  of  temperature, 
barometric  pressure,  relative  humidity,  etc. 

2)  If  any  two  objects  are  put  on  the  pan,  the  indicator  reading  is  the  sum  of 
the  readings  for  each  object  placed  individually. 

The  first  property  means  that  the  imprecision  of  the  balance  readings  is  zero.  The 
second  property  means  that  the  balance  is  perfectly  linear.  Later  we  will  show  how 
to  use  real  balances  which  only  approximate  these  two  important  features. 

A  final  property  of  our  idealized  balance,  which  it  shares  with  actual  scales 
and  balances  is  that  the  indicator  responds  to  a  force  on  the  pan,  not  a  mass.  In 
other  words,  our  balance  would  give  different  readings  on  the  moon  than  on  earth 
(though  properties  1)  and  2)  would  be  unaffected). 

When  the  balance  has  reached  an  equilibrium  condition,  the  sum  of  all  forces 
acting  on  the  system  must  be  zero.    We  then  have 


and  since  the  forces  acting  on  the  system  are  known  to  act  in  one  direction  and  its 
exact  opposite  (i.e.  up  and  down)  we  can  once  again  ignore  the  vector  notation  and 
wri  te 


n 


i=i  1 


n 


Ef-o 

i=l  1 


(4) 


where  n  =  total  number  of  forces  acting  on  the  system. 


The  forces  acting  on  the  object  placed  on  the  pan  of  the  balance  are: 


Ft  =  Mg 


gravitational  force 


F2,    balance  force  exerted  by  the  balance  via  deflection  of  a  beam,  stretching  of  a 
spring,  or  some  other  method.    Obviously  this  force  has  to  be  opposite  to  F, . 
Since      acts  opposite  to  F-j ,  we  adopt  the  convention 
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Figure  1.    An  idealized  scale.    The  indicating  dial  is  marked  in  equally  spaced 
increments.    The  number  pointed  to  by  the  needle  can  be  resolved  to  any  desired 
precision. 


where  6,  is  the  balance  reading  and  k  is  a  positive  constant  (yet  to  be  evaluated) 
which  converts  balance  readings  to  mass  units. 

In  addition  to  these  main  forces,  we  have  to  consider  buoyant,  electrostatic, 
and  electromagnetic  forces  on  the  object.    Of  these  three  only  the  first  one  plays  a 
significant  role  in  usual  metrology  and  we  now  look  at  the  effect. 

Under  the  influence  of  gravity,  the  force  on  a  body  submerged  in  a  fluid  is 
equal  in  magnitude  but  opposite  in  direction  to  the  weight  of  the  displaced  medium. 
This  upward  force  is  known  as  the  "buoyant  force."    From  the  definition  we  see  that 

F3  =  -mg 

where  m  =  mass  of  displaced  medium.    The  buoyant  force  was  first  quantified  by 
Archimedes  in  the  third  century  B.C.    We  also  can  write 

m  =  Vpm 

,  Km 
so  that 

F3  =  -Vpmg  (6) 

where 

V  =  volume  of  the  displaced  medium 

=  volume  of  object 
Pm  =  density  of  medium  =  m/V  . 

With  these  three  forces,  we  can  now  write  eq.  (4)  as 

Fl  +  h  +  F3  ■  0 

or 

F1+F3=-F2  . 
Substituting  from  eqs.  (2),  (5)  and  (6)  yields 

Mg  -  pmv-g  =  ke^ 

(M  -  PmV)g  =  ke^  (7a) 

(M  -  PfnV)  =  ke1    .  (7b) 

It  is  very  important  to  realize  from  eq.  (7a)  that  if  the  weighing  is  done  in  a 
vacuum  chamber  then  the  F3  term  drops  out  and  eq.  (7a)  reduces  to 

Mg  =  ke^    .  (7c) 

Furthermore,  the  medium  or  environment  can  be  any  fluid--air,  water,  oil,  etc. 
For  most  applications  and  for  practical  purposes,  the  medium  is  usually  air.  Even 
though  air  is  a  mixture  of  gases  and  the  density  is  quite  small  compared  to  densities 
of  objects  usually  measured,  the  correction  term,  F3,  cannot  be  ignored  by  metrolo- 
gists,  as  we  shall  see. 

Let  us  now  consider  the  weighing  of  two  objects  of  mass  Ml  and  M2  under  the  same 
conditions,  i.e.  in  the  same  medium.    We  then  have  from  eq.  (7a) 


or 


g(M1-PmV1)  =  ke^  (7d) 
g(M2-pmV2)  =  k92g    .  (7e) 
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The  density  of  an  object  is  its  mass  devided  by  its  volume.    That  is, 

where  p-|  is  the  density  of  object  1.    We  can,  therefore,  always  replace  a  quantity 
such  as 

Ml  -  eJl 

by  the  identical  quantity 

M,(l  -  Pm/Pl)  . 

The  choice  of  representation  is  completely  optional,  although  the  proper  choice  can 
often  simplify  a  calculation.    For  the  present,  we  will  continue  with  the  volume 
representation. 

Subtracting  eq.  (7e)  from  7(d)  yields 

9(M1  -  p^)  -  g(M2  -  pmV2)  =  g  k6-|  -  g  ke2  (7f) 

and  by  setting 

e,  -  e2  =  A6 

we  have 

9(M1  "  pmV  "  9(M2  "  pmV2)  =  9kAe    •  (8) 
Cancelling  g  yields 

M1  -  M2  =  kA6  +  pm  (V1  -  V2)    .  (9) 

Equation  (9)  is  one  of  the  fundamental  equations  in  mass  metrology  and  is  known 
as  the  "Weighing  Equation."    Its  importance  will  become  evident  in  Section  3.  For 
now  we  simply  point  out  that,  using  the  best  balances  available,  one  can  determine 
A6  to  much  higher  accuracy  than  either  6-j  or  02-    Thus  eq.  (9)  turns  out  to  be  more 
useful  than  eq.  (7b) . 

In  most  cases,  the  laboratory  conditions  are  such  that  the  medium  is  air,  in 
which  case 

PA  =  Pm  =  density  of  air 

and  eq.  (9)  becomes 

M1  -  M2  =  kA6  +  pA(V1  -  V2)  (10) 

Note:        When  we  refer  to  "weighing"  in  this  text,  we  will  always  mean  the  process 
which  determines  the  product  k  •  6.    It  should  be  clear  from  the  above 
paragraphs  that  finding  the  mass  of  an  object  involves  more  than  just 
weighing  the  object.    Also,  we  have  no  future  need  to  refer  to  "weight"  as 
defined  by  (3).    From  this  point  on,  the  term  "weight"  will  be  reserved 
solely  to  designate  an  object  manufactured  to  have  particular  nominal 
mass.    For  example,  a  "1-g  weight"  is  an  object  whose  mass  is  close  to  1  g. 
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From  eqs.  (9)  and  (10)  it  is  obvious  the  condition 


can  exist  if  and  only  if 


PJV,  -  v2)  e  0  . 

This  can  only  hold  if 

1 )  p    =  0  i.e.,  the  measurement  is  done  in  vacuum 


or 


2)  V,  =  V2 


The  first  case  is  usually  not  encountered  in  metrology.    The  second  case  can  be 
approximated  quite  easily.    Consider  for  example  two  weights  of  nearly  equal  mass, 
such  that 

M1  ~  M2  , 

made  from  the  same  billet  of  alloy.    In  this  case,  the  densities  are  equal  so  that 

Pi  =  P2  =  p  • 


Since 


M,  M9 
V1  =  -J-  and   V2  =  — 

pl  p2 


we  have  from  (10) 


n 

."-1 

\pl 

p2, 

-H2 

M1  -  M2  =  kA6  +  pA 

(M1  -  M2)  (1  -  — )  =  kA6 

P 

M        m  kA6 
M-J  -  M2  =   


l-PA/p 


M1  -  M2  =  kA6 


because  pA  «  p  and  M-j  -  M2  is  already  known  to  be  a  small  number. 

Relations  such  as  that  expressed  by  eq.  (1)  or  (10)  are  of  use  to  metrology 
only  if  all  variables  are  expressed  in  a  consistent  set  of  units.    In  what  follows 

we  will  assume  that  we  are  working  in  the  International  System  of  Units 
(also  known  as  the  Systeme  International  des  Unites,  or  simply,  the  S.I.). 
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The  unit  of  mass  in  the  S.I.  is  the  kilogram.    By  definition,  one  kilogram  is 
exactly  equal  to  the  mass  of  an  object  known  as  the  International  Prototype  Kilogram, 
or  I.P.K.    This  object  is  made  of  a  platinum-iridium  alloy  and  is  stored  at  the 
International  Bureau  of  Weights  and  Measures  in  Sevres,  France.    All  laboratory  mass 
standards  must  ultimately  be  traceable  to  the  I.P.K. 

2.1.    True  Mass  and  Apparent  Mass 

So  far,  the  "mass"  we  have  considered  is  that  found  in  Newton's  basic  equation 
(i.e.,  eq.  (1)).    For  that  reason,  this  concept  of  mass  is  sometimes  called  "true 
mass".    It  is  a  quantity  intrinsic  to  the  object  alone  and  not  to  its  situation.  We 
have  also  seen  that,  in  a  vacuum,  the  force  registered  by  the  balance  is  exactly 
proportional  to  the  Newtonian  mass.    For  this  reason  the  "true"  mass  is  also  known  as 
the  "vacuum"  mass.     The  terms  "mass",  "true  mass"  and  "vacuum  mass"  are  interchange- 
able. 


The  need  to  apply  buoyancy  corrections  to  mass  measurements  has  led  to  the 
adoption  of  so-called  "apparent"  masses.    As  explained  below  and  in  Appendix  A,  use 
of  apparent  mass  can  simplify  the  job  of  both  the  weight  manufacturer  and  the  metro- 
logist.    The  drawback,  however,  is  that  one  must  learn  a  new  concept  (that  is, 
apparent  mass) . 

Of  limited  scientific  significance,  the  apparent  mass  approach  nevertheless 
often  simplified  calculations  for  the  metrologist  at  a  time  when  exact  computations 
were  tedious  (either  done  by  hand  or  through  nomographs)  and  difficult  to  check. 
Today,  when  every  metrologist  has  access,  at  the  very  least,  to  a  digital  pocket- 
calculator,  the  argument  for  apparent  masses  is  less  compelling.    Nevertheless,  both 
for  completeness  and  because  commercial  weights  are  still  manufactured  on  the  basis 
of  apparent  mass,  we  cannot  avoid  this  topic. 

The  fact  is  that  most  mass  metrology  is  carried  out:    (1)  in  air;    (2)  in 
laboratories  near  sea  level;  and  (3)  at  a  temperature  of  ^20  °C. 

3 

Air  at  such  conditions  has  a  density  nearly  equal  to  1.2  mg/cm  .  Fluctuations 
about  the  mean  density  are  known  rarely  to  exceed  3  percent.    Thus  97  percent  of  the 
buoyancy  correction  can  be  "built  into"  the  calibration  of  a  weight  so  that  small 
deviations  from  the  correction  are  all  that  need  be  considered. 

Let  us  carry  through  an  example  by  defining  a  reference  metal  of  density  pn. 

the  true  mass  of  an  object  made  of  reference  material  R. 
(This  object  need  not  actually  exist.) 


Let 


K  = 


pR  =  the  density  of  the  reference  material 

1.2  mg/cm3  =  the 
nearly  the  densi' 
mm  Hg  pressure. ) 


pn  =  1.2  mg/cm   =  the  reference  density  of  air.    (This  is  very 
nearly  the  density  of  air  at  20  °C,  50%  humidity,  and  760 


tn  =  20  °C  =  reference  temperature. 


If  one  actually  were  to  weigh  an  object  under  vacuum  conditions  to  determine  its 
mass,  one  would  have  to  be  certain  that  the  object  weighed  was  sufficiently  stable- 
that  is,  no  gasses  normally  adsorbed  on  the  surface  at  atmospheric  pressure  are 
vacuumed  off,  the  weight  itself  does  not  "out-gas,"  etc.    The  equality  between  vacuum 
mass  and  true  mass  assumes  that  the  object  being  weighed  is  stable. 
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We  now  define  the  "apparent"  mass  of  an  object  as  follows:    The  "apparent"  mass 

of  an  object  X  is  equal  to  the  "true"  mass  MR  of  just  enough  reference  material  to 

produce  a  balance  reading  equal  to  that  produced  by  X  if  the  measurements  are  done  at 
temperature  tg  in  air  of  density  pQ  (see  fig.  2). 

Stated  another  way,  we  first  define  specific,  unique  weighing  conditions: 

3 

i.     the  air  density  is  1.2  mg/cm 

ii.     the  temperature  is  20  °C  (it  is  necessary  to  specify  temperature  because 
the  volume  of  a  weight  depends  slightly  on  temperature.) 

Next  we  imagine  a  mass  comparison  between  an  object  X  and  an  assembly  of  known 
objects  R.    The  definition  specifies  that  the  R  objects  shall  all  have  a  density 
equal  to  pp.    The  mass  of  the  R  objects  is  then  adjusted  until 

A6  =  eR  -  ex  =  o  . 

When  these  conditions  are  all  satisfied,  MR  is  by  definition  equal  to  the  apparent 

mass  of  X.  We  never  really  have  to  carry  out  the  experiment  because  from  eq.  (9)  we 
have 

MX  -  MR  =  kA6  +  p0  (VX  "  V 

where  VR  is  given  by  VR  =  Mr/p^-    We  can  rewrite  this  relation  using  density  nota- 
ti  on : 


so  that 


and  with 


MX  "  PQVX  "  MR  +  p0VR  =  kA6 
(1  -  P0/PX)  -  Mj  (1  -  p0/pr)  = 


kA0 


A6  =  0 


T         My    (1-Pn/Py)  « 

Ml     =    -2       =     My  .  (11) 

(1-Pq/Pr) 


Note  that  the  denominator  is  simply  a  constant  whose  value  is  the  same  for  any 

3  A 
weight,  X.    Also,  for  px  greater  than  Ig/cm    (the  density  of  water),      will  be 

within  0.2  percent  of  mJ.  Furthermore,  eq.  (11)  states  that  the  apparent  mass  can  be 
obtained  by  multiplying  the  true  mass  by  a  buoyancy  factor. 

The  usefulness  of  apparent  mass  and  hence  eq.  (11),  as  well  as  the  subtle 
distinctions  between  true  and  apparent  mass  are  shown  in  the  next  five  subsections. 
We  begin  with  an  example. 


9 


QC 
00 

h- 
x 
g 

LJJ 


t—  Of 

II 

00 

00 

MA 

cc 
Q_ 

LU 

■  I 

II 

ZD 

>- 

DC 

1- 

00 

_J 

z 

< 

LU 

1- 

Q 

O 

y— 

I 

o 

LU 


oo 

00 

< 


LU 

ZD 

cc 


Q_ 

HO 
> 

tz 

oo 
z 

IULJI 


> 


QC 

IUN 

II 

Z 

LL 

<D 

LU 
Q 

>■ 

m 

Z 

z 

LU 

m 

I 

x 

l!  - 

oo 

00 


LU 
ZD 
DC 

t— ■ 

OD 

X 


00 
00 

< 


z 

LU 

< 

□l 

QL 
< 


QJ 

JZ 

+J 

C 

•  1 — 

(/) 

03 

E 

-M 

c 

ai 

s_ 

O 

03 

•r— ■ 

a.  +j 

Q. 

03 

C 

■  1 — 

f- 

<+- 

O 

QJ 

T3 

C 

O 

QJ 

jz 

+J 

+-> 

•  1 — 

sz 

c 

•1 — 

4- 

OJ 

TD 

~o 

QJ 

•1 — 

QJ 

4- 

JZ 

'1  

+J 

(_> 

QJ 

4— 

Q- 

O 

CO 

C 

CO 

O 

•  1 — 

+-> 

SZ 

03 

O 

+J 

C 

4-> 

ai 

03 

c/i 

E 

<u 

S_ 

i. 

O 

Q.  4— 

QJ 

sz 

S- 

■1 — 

u 

- 1 — 

aj 

+J 

X 

ro 

0 

E 

xj 

QJ 

_SZ 

QJ 

O 

JZ 

CO 

h- 

< 

CO 

CO 

03 

QJ 

5- 

Z3 

CC 

CT) 

-J. 

Ll_ 

10 


2.1.1.  Example 


Using  the  ideal  balance  described  earlier,  we  wish  to  measure  the  mass  of^an 
unknown  weight,  X,  in  terms  of  standard  weights,  S,  of  known  apparent  masses  M$'s. 
During  the  actual  measurement,  done  in  air  of  temperature  t  and  air  density  p^,  it 
was  found  that  a  summation  of  standards  produced  the  same  reading  on  the  balance  as 
the  unknown  mass. 

Under  what  conditions  can  we  set 

My    «  ? 
A  O 

Sol ution: 

Starting  with  the  balance  equation  for  two  masses  (eq.  (10)). 

M]  -  M2  =  kA6  +  pA  (V1  -  V2) 

we  have 

mJ  -  h\  -  k(ex  -  es)  +  pA(vx  -  vs) 
■  pa<vx  "  V 

since  the  measurement  yielded  6X  =  6<.. 
Finally,  we  can  write  this  as 

MX  "  MS  =  PAAV    '  (E1) 
Similarly,  we  have  for  the  apparent  masses,  using  eq.  (11): 


By  setting 


we  simplify  to 


(1-p0/pr)  (l-p0/pR) 

T  T  Y  ^ 

K  -  MS»  -  j  PQ  +  pf  p0 
(1-P0/PR) 


-'ho  •  -=vso  at20°c 

px  ps 


HA,  HA=K-^-PQ  'VXO-W     .  (£2) 

X  0-p0/pr) 
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Substituting  for      -       from  (El),  we  have 


MA     MA     PA(VVS}  "  p0  (VX0"VS0) 
(1-Pq/Pr) 

Volumes  have  a  slight  temperature  dependence  but  if  the  temperature  of  the 
measurement  is  close  to  20  °C  (generally  true  in  the  laboratory),  then 


and 


vx  *  vxo 


V    «  V 

vs  vso 


so  that  (E3)  reduces  to 

(PA-p0)(W=     ApAV  .  (E4) 

(1-Pq/Pr)  (1-Pq/pr) 

The  term  1/(1-Pq/pr)  can  be  approximated  by  1  for  the  following  reason:    Via  a 
Taylor  expansion 


1-p0/pr        pr    Vpr/  Vp 


Now  since  pQ  is  the  density  of  air  and  pR  is  the  density  of  a  reference  metal,  we 
have 

p0  -3 
—  <  10  J 


so  that 


R 


-  =  ApAV  +  ApAV  (p0/pR)  +  . . • 


(1-P0/PR) 


and  since  ApAV  is  already  a  small  number,  a  correction  of  10     (or  0.1  percent)  is 
insignificant.    Hence  only  the  first  term  has  to  be  retained.    This  argument  is 
important  in  that  it  is  used  time  and  again  in  estimating  buoyancy  corrections. 

With  these  approximations,  we  finally  arrive  at  the  equation  for  apparent  mass 
differences 

My  -       =  ApAV    .  (E5) 
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We  now  would  like  to  find  the  conditions  for  which 


MX*MS 


This  implies  that  ApAV  «  0. 


When  one  encounters  such  a  relation,  one  must  always  decide  its  precise  meaning. 

In  this  particular  case  the  meaning  is:    ApAV  is  smaller  than  the  acceptable  un- 
A 

certainty  of  M„. 

A 

For  example,  suppose  Mx  needs  to  be  known  to  one  part  per  million.  Then 


really  means  that 


ApAV  ~  0 


ApAV 


M 


<  10 


-6 


(E6) 


A  T 
But  we  already  know  that      does  not  differ  from      by  more  than  a  fraction  of  a 

percent  (see  comments  after  eq.  (11)).    We  can  therefore  approximate 


X 

Pv  =  ~ 


by 


Now  (E6)  becomes 


Pv  * 


Ap 

AV 

<  10-6 

pX 

VX 

(E7) 


The  volume  ratio  AV/V^  can  be  written  as 

AV  _  Yvs 

VX  V 

In  the  case  of  different  alloys  of  stainless  steel  (SS),  VY  will  be  close  to  V<»,  so 
that 


AV 


<  10 


-1 
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We  now  have 


x  TO"1  <  10"6 


(E7') 


The  density  ratio  can  be  written  as 


Ap  .  PA"P0 
PX  PX 

3 

p^  *  8.0  g/cm    for  SS 


where 
and 

so  that 


1.2  x  10"3  g/cm3 


Ap_  (pA-1.2xlO-) 
PX  8.0 

The  inequality  (E7')  now  reduces  to 

^-3 


PA' 


1.2x10 


8.0 


x  10"1  <  10"6 


which  is  satisfied  as  long  as 


1.28  x  10"3  g/cm3  >  pA  >  1.12  x  10"3  g/cm3. 

This  indicates  that  the  air  density  during  the  measurement  can  be  off  by  as  much  as 

3 

0.08  from  1.20  mg/cm    and  still  satisfy  condition  (E6). 

In  contrast,  the  buoyancy  correction  for  true  masses  is  much  more  critical. 
From  eq.  (El),  we  want  to  satisfy  the  condition  that 


PAAV 


<  10 


-6 


or 


A 


AV 


<  10 


-6 


Using  the  same  material  as  before,  we  have 


so  that 


—  x  10"1  <  10"6 
8.0 


PA  <  0.08  x  10'3  g/cm3 


14 


Note  that  in  order  to  keep  the  desired  accuracy,  the  air  density  during  the  measure- 
ment can  only  be  ignored  if  it  is  very  small,  i.e.  vacuum. 


This  has  been  a  long  example  in  which  several  generally  useful  techniques  have 
been  introduced.    It  is,  therefore,  worthwhile  to  look  back  and  summarize  what  this 
example  also  teaches  about  apparent  mass.    The  most  important  inference  we  can  draw 
is  that,  for  true  mass  comparisons,  we  almost  always  need  to  account  for  air  buoyancy 
even  though  we  might  be  able  to  use  an  assumed  value  for  the  air  density  (recall  that 
air  density  is  usually  constant  in  a  given  laboratory  to  about  ±3  percent).  For 
apparent  mass  comparisons,  the  air  buoyancy  correction  is  much  less  important  and 
can,  therefore,  sometimes  be  ignored.    Apparent  mass  works  the  way  it  does^because 
the  formalism  already  has  an  assumed  value  for  air  density  (i.e.  1.2  x  10~w  g/cm^) 
built  into  it.    In  measurement  situations  where  this  assumed  value  is  not  accurate 
enough,  the  buoyancy  term  in  eq.  (E5)  becomes  important.    This  eventuality  occurs  at 
precisely  the  same  point  when  pfl  in  eq.  (El)  can  no  longer  be  sufficiently  well 

-3  3 

approximated  by  1.2  x  10        g/cm  . 

2.1.2.    Reference  Materials 

So  far  we  have  not  yet  defined  a  reference  material.    At  present,  two  different 
apparent  mass  bases  are  utilized  by  the  NBS.    One  of  them,  the  older  one,  is  called 
"normal  brass"  and  was  the  logical  choice  when  most  laboratory  weights  were  made  of 
brass. 

Normal  brass  is  defined  by 

PBQ  =  8.4  g/cm3  at  0  °C 
-5 

aB  =  5.4  x  10   /°C  =  coefficient  of  cubical  expansion 
3 

p«  =  1.2  mg/cm   =  air  density  at  20  °C. 

These  parameters  together  with  the  expression  for  volumetric  expansion 

VB(t)  =  VBO[l+aB(t-0  °C)]  (12) 

and  the  corresponding  density  correction 

PB(t)  =  PBO/[l+aB(t-0  °C)]  (13) 

completely  determine  the  parameters  at  20  °C.    Note  from  above  that  pBq  is  given  at 
0  °C  and  therefore  has  to  be  determined  for  20  °C  via  eq.  (13)  to  yield 

pR(20)  =  =  8.390938  g/cm3    .  (13') 

b  [!+5.4xl0"bx20] 

The  second  apparent  mass  basis  is  referred  to  an  arbitrary  material  with  the 
following  density: 

pA  =  8.0  g/cm3  at  20  °C  . 

Once  more, 

p0  =  1 .2  mg/cm3  at  20  °C. 

Note  that  this  basis  does  not  require  any  expansion  coefficient  since  all 
parameters  are^  defined  at  20  °C.    It  is  therefore  comparable  to  the  basis  defined  for 
normal  brass  via  eq.  03'). 
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It  has  now  become  common  practice  to  report  masses  and  their  uncertainty  in  the 

true  mass  basis  and  then  quote  corrections  to  nominal  values  for  "normal"  brass  and 

3 

for  pA  =  8.0  g/cm    bases.    Clearly,  the  latter  apparent  mass  basis  is  beginning  to 
be  preferred  by  many  laboratories  since  most  weights  used  in  a  set  are  made  of 
stainless  steel  which  has  a  density  very  close  to  8.0  at  20  °C.    Thus  the  correction 
term  in  eq.  (11)  is  small,  making  the  apparent  mass  close  to  the  true  mass  (see 
Appendix  A. ) 

We  are  now  in  a  position  to  write  eq.  (11)  in  the  "normal  brass"  (pb)  and 
in  the  "8.0"  (pg)  basis.    Starting  with  eq,  (11),  we  have 

MA     MT  (1-p0/pX} 
O-Pq/Pr) 

and  substituting  for  p^  with  p^  from  (13),  we  obtain 


MX  =  MX 


P0[l+aB(t-0  °C)] 
(1-  T 
PB0 


or 

ha  .  mt       0-P0/PX>  (14) 

pn 

1  -  [l+20aB] 
PB0 

where  the  denominator  is  just  a  constant  equal  to  0.9998569886. 
Similarly,  we  use 


3 

for  the  8.0  g/cm    basis  to  yield 


PR  =  p0 


A     __    T  (1-Pp/Px) 


X       X  (1-p0/p6> 

where  the  denominator  is  once  again  a  constant  exactly  equal  to  0.99985000.    It  is 
interesting  to  note  that  apparent  masses  in  the  two  bases  we  have  examined,  eqs. 
(14)  and  (15),  have  a  constant  ratio 

(M  )ddacc 

 =  1.00000699  . 


(MA) 


8.0 


That  is,  the  apparent  masses  of  the  same  weight  in  the  two  bases  differ  by  only 
^.0007  percent. 
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Equation  (15)  indicates  a  verv  important  point.  If  the  object  being  weighed 
has  a  density  px  equal  to  8.0  g/cm    (i.e.  px  =  pi)  at  20  °C  then 


For  the  brass  basis,  eq.  (14),  a  similar  argument  has  to  be  treated  with  more 
caution.    In  this  case,  if  px  is  equal  to  the  density  defined  by  eq.  (131),  then 


M  = 


2.1.3.    Reporting  Apparent  Mass 


Sometimes  the  apparent  mass  correction  in  a  given  basis  is  required.    This  is 
defined  as 

Correction  (mJ)  =      -  Nx  (16) 

where       is  the  nominal  value  of  the  weight. 

2.1.4.    Comparing  Two  Apparent  Masses 

We  now  derive  an  equation  similar  to  (10),  appropriate  for  apparent  mass  differ 
ences.    Using  eq.  (11)  we  have  for  mass  1 

MA  =  MlT-p0V01 


1-Pq/Pr 


where  we  have  used  mT/p-|  =  Vq-j  . 
Similarly  we  have  for  mass  2 


MA  M2T-p0V02 
M2  =   


1- 


P0/PR 


Subtracting  the  second  from  the  first  yields 

1  ]-Po/pR 
From  eq.  (10),  we  have  the  true  mass  difference  of  the  two  weights 

M ]  -      =  kA6  +  pA  (VrV2) 

so  that,  upon  substitution,  we  obtain 

MA  „A  kA9+pA(V1-V2)-p0(V01-V02)  _  (17) 
Mq  -  M«  =  ■  

Now,  employing  eq.  (12),  we  can  use 
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Vl  ■  V01  H+a^t-to)] 

and 

V2  =  VQ2  [Ha2(t-t0)J 
where  Vq-j  and  Vq2  are  the  respective  volumes  at  tg  =  20  °C. 
Substituting  for  V-j  and  V2  into  (17)  yields 

-       =  kAO  -  p0  (VQ1-V02)  +  pA  [V01+V0iai(t-20)  -  VQ2-V02a2(t-20)] 


where  we  have  set 


tQ  =  20  °C 


and  have  approximated  (1-Pq/Pr)"^  by  1  (see  derivation  of  (E5)). 
Simplifying  further  gives 

-  ^  =  kA6  +  (pA-P0)(v0rV02)+pAC(t"20)(V01arV02a2)]     ■  (18) 

Equation  (18)  is  a  fully  corrected  formula  for  calculating  mass  differences  in  the 
apparent  mass  basis,  except  for  the  one  approximation  mentioned  above.     At  this 
point  it  is  worthwhile  to  consider  a  specific  example  in  order  to  get  a  feeling  for 
the  magnitude  of  the  terms  in  eq.  (18).    Assume  that  two  weights  of  nominal  value  2g 
have  the  following  properties: 

VQ1  =  0.2564  cm3 

VQ2  =  0.1884  cm3 

also,  a1  a  o2  =  4.5  x  10"5/°C; 

3 

PA  =  1.17  mg/cm 


3 

Pq  =  1.2  mg/cm  , 


We  then  find  that 


(PA  "  P0)(V01  -  V02)    =  -0.03  x  0.068  =  -2,04  x  10"3  mg 


and 


PA(t  -  20)(V01a]  -  VQ2a2)    =  117  (t  -  20)  x  3.06  x  10"6 

=  (3.58  x  10"6)(t  -  20)  . 
Considering  an  extreme  condition  such  that  t  is  10  °C,  we  have  for  the  latter  term 

=  -3.58  x  10'5  mg  . 
Combining  this  term  with  the  first  we  have 

-2.04  x  10"3  -  0.0358  x  10"3  =  -2.08  x  10"3  mg 
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so  that  the  second  term  (under  extreme  conditions)  contributes  only  0.04  yg  to  the 
overall  buoyancy  correction.    Usually  the  (t  -  20)  term  is  of  the  order  of  0  °C  to 
4  °C  so  that  the  effect  of  the  second  term  on  the  overall  buoyancy  correction  of  the 
above  weights  is  even  less.    Hence  in  most  cases  the  second  term  can  be  safely 
ignored  and  eq.  (18)  can  be  approximated  by 

-       =  kA6  +  (pA  -  p0)(V01  -  VQ2)  .  (19a) 

For  true  masses  we  start  with  eq.  (10)  and  write 

m{  -       =  kA6  +  pA  (VrV2)  . 

Employing  the  same  volume  expansion  formulae  used  for  the  apparent  masses  (after  eq. 
(17) ) ,  we  derive 

m{  -  M2  =  kA6  +  pA  (VQ1  -  V02)  +  PA(t-tQ) 

x  (VQ1a1  -  VQ2a2)    .  (10') 

Once  again  the  second  term  is  very  small  compared  to  the  first  (see  arguments  above) 
and  the  relation  can  usually  be  approximated  by 

M{  -  M2  =  kA6  +  pA  (VQ1  -  V02)    .  (19b) 

A  quick  comparison  between  eqs.  (19a)  and  (19b)  shows  that  the  apparent  mass  differ- 
ence has  a  correction  term  which  contains  the  factor  (pA  -  Pq)  whereas  for  true 
masses  the  corresponding  factor  is  p/\.    This  then  shows  that  the  former  correction 
is  much  smaller.    In  both  cases  the  unknown  parameter  is  pA- 

2.1.5.    Sources  of  Error  in  Buoyancy  Corrections 

It  is  important  to  determine  the  sources  of  error  in  buoyancy  corrections  in 
order  to  make  the  best  measurements  possible  with  the  least  amount  of  effort  and 
expense.    Let  us  look  back  to  eqs.  (19a)  and  (19b)  and  write  the  most  important 
buoyancy  terms: 

True  Mass  Apparent  Mass 

pA(V01  "  V02}  (pA  -  P0)(V01  "  V02} 

Since  pq  is  defined  exactly,  uncertainties  in  the  mass  difference  measurements 
arising  from  buoyancy  corrections  are  due  to  errors  in  pa,  Vq-j  ,  and  Vq2-    Let  us 

denote  the  uncertainty  in  a  quantity  by  a  preceding  6  (i.e.,  the  uncertainty  in  pA 

is  6pA).  Then,  from  eqs.  (19a)  and  (19b)  the  sources  of  uncertainty  in  the  buoyancy 
corrections  are 

True  Mass  Apparent  Mass 

6pA(VQ1  -  VQ2)  6pA(VQ1  -  VQ2) 

6VQ1  pA  6VQ1(pA  -  pQ) 

6VQ2  pA  6VQ2(pA  -  p0) 

where  we  will  assume  that  6VQ1 ,  6VQ2,  and  <5pA  are  uncorrelated. 
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In  the  calculation  of  true  mass  differences,  uncertainty  in  the  volume  of  the 
masses  is  much  more  important  than  it  is  for  differences  in  apparent  mass  as  can  be 
seen  by  comparing  similar  terms: 

(6V01  ♦  6V02)pA  >  («¥01  +  6V02)(pA-p0)  . 

Also  we  see  that  any  uncertainty  in  Vq-j  or       will  lead  to  systematic  uncertainty 
in  mass  measurements  (even  if  Vq-j  =  VQ2). 

Example: 

Two  weights  produce  the  same  reading  on  an  ideal  balance  in  air  of  p^  =  1.17  mg/cm 
and  t  =  20  °C.    Both  weights  are  said  to  be  of  density  7.89  g/cm3.  Find 

1)  their  true  mass  difference 

2)  their  apparent  mass  difference  against  Pq  =  8.0  g/cm 

3)  the  errors  in  these  numbers  due  to  uncertainty  in  the  volumes. 

The  weights  are  marked  "1  kg". 
Answers: 

1)  From  eq.  (19b) 

t  Pa  T  pA 

M      (1-  -A- )  -  M'     (1-  =  k-A6  =  0 

1  7.89        c  7.89 

because  A6  is  given  as  0. 
Hence  m|  -  mJ  =  0. 

2)  From  eq.  (19a) 

A  "  M2  =  (pA  "  po)(V01  '  V02}  +  kAe  =  0 

since  A6  =  0  and        =  V^^  (because  mJ  =      and  p^  =  P^)- 

3)  the  weights  could  have  come  from  different  lots  of  metal.    Since  no  uncer- 
tainty is  given  for  density,  we  must  assume  that 

=  7.89  ±  0.005  g/cm3  , 

i.e.  the  density  is  uncertain  in  the  first  unreported  decimal  place. 

.As  a  worst  case,  let 

Pq-j  =  7.895  and  p^  =  7.885  . 

Then 
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6(M 


,T 


7.885 


Ml 


7.895 


7.885 


7.895 


) 


0.00117  x  1000 


a  0.270  mg 


5(W 


■A 


0.00003  x  1000 


-  0.005  mg 


Note  that  the  nominal  mass  value  of  each  weight  is  accurate  enough  to  use  in  esti- 
mating the  uncertainties  in  question. 

CAUTION:    In  calibration  certificates,  it  is  normal  practice  not  to  include  in 
the  assignment  of  systematic  errors  any  contribution  from  volume  uncertainties:  The 
volumes,  usually  calculated  from  densities  provided  by  the  weight  manufacturers,  are 
assumed  to  be  without  error.    For  the  most  critical  work,  volumes  of  individual 
weights  must  actually  be  measured.    (Since  1982,  NBS  has  been  routinely  measuring  the 
volume  of  single-piece  1 -kilogram  weights  which  require  the  best  possible  mass 
calibration.) 

We  have  yet  to  examine  the  effects  on  mass  measurements  of  uncertainty  in  p.. 
For  mass  differences,  we  have  shown  above  that  the  errors  are  approximately 
6p^(Vq.|  -  V-Qg).    Therefore,  unless        =  Vq2»  we  mst  be  concerned  with  uncertainty 

in  p^.    The  calculation  of  p^  along  with  an  analysis  of  the  attendant  uncertainties 

are  presented  in  the  next  section. 

2.2.    Air  Density 

We  have  seen  from  the  derivations  in  Section  2.1  that  the  correction  to  the  mass 
determination  requires  the  precise  knowledge  of  the  air  density  (p^)  if  the  weighing 

is  done  in  an  air  environment.  (We  will  assume  from  now  on  that  we  are  always  work- 
ing in  air.)    The  buoyancy  correction  in  eq.  (19b)  requires  p^  for  true  masses, 

whereas  for  apparent  masses  we  require  (p^-Pg)  in  ecJ«  (19a).    It  should  be  remembered 

that  Pq  is  defined  as  pQ  =  1.2  mg/cm^  which  is  the  approximate  value  for  20  °C  with 

1  atmosphere  of  pressure  with  50%  relative  humidity. 

The  usual  approach  to  determining  p.  is  to  calculate  the  density  of  air  based  on 
an  "equation-of-state."    To  use  this  equation,  one  must  supply  ambient  values  for 
temperature,  barometric  pressure,  and  relative  humidity.    For  extremely  exacting 
work,  the  concentration  of  carbon  dioxide  in  the  laboratory  ambient  is  also  required. 
Appendix  B  contains  a  derivation  of  the  equation-of-state  based  on  the  work  of  Jones 
[2J.    This  formulation  is  close  to  those  which  are  found  in  handbooks  or  older 
references  but  has  three  unique  virtues: 
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1.  It  is  based  on  the  most  recent  reference  data  available. 

2.  It  has  been  derived  in  such  a  way  that  uncertainties  associated  with  the 
equation  itself  are  known  and  stated. 

3.  Virtually  the  same  equation  has  been  endorsed  for  international  use  by  the 
Consultative  Committee  for  Mass  of  the  International  Committee  of  Weights 
and  Measures  (C.I. P.M.). 3 

The  interested  reader  is  urged  to  consult  Appendix  B  and  especially  Ref.  [2]  for 
the  complete  derivation.    We  state  here  that4 

(ms_)  =  0-0034848    .  (p  _  0>0037960  .  u  .  eJ  (20a) 
A       3       (t+273.15)  s 


where        t  =  temperature  in  °C 

P  =  barometric  pressure  in  pascals 

(133.322  Pa  =  1  mm  Hg) 
U  =  relative  humidity  in  % 
e^  =  saturation  vapor  pressure. 

e$  =  1.7526  x  1011  exp  (-5315.56/(273.15  +  t))Pa 
is  an  excellent  approximation. 


If  barometric  pressure  is  measured  in  mm  Hg,  then 


(!02_)  =     0- 46460    (P.0.0037960  u  .  e j  (20b) 
3       (t+273.15)  s 


cm 


o 

where  es  =  1.3146  x  10  exp  (-5315.56/(273.15  +  t))  mm  Hg  is  an  excellent  approxi- 
mation. 


2.3.    Sensitivity  Arguments 

The  major  correction  normally  required  for  mass  determination  is  due  to  air 
buoyancy.    As  is  evident  from  eqs.  (19a)  and  (19b)  the  only  parameter  that  has  to  be 
determined  very  carefully,  aside  from  the  volumes  of  the  weights  being  compared,  is 
the  air  density.    Equation  (20)  clearly  indicates  that  three  parameters  affect  the 
air  density,  namely  humidity,  barometric  pressure,  and  temperature.  We  will  see  in 
Section  2.3.5  that  an  overall  uncertainty  in  p«  of  ^0.0030  mg/cm  ,  i.e.  ^0.25  percent 

of  p^  near  sea  level,  is  usually  acceptable.    In  Appendix  C  we  see  that  the  total 

uncertainty  in      is  the  square  root  of  the  sum  of  the  squares  of  the  individual 

uncertainties  arising  from  measurements  of  humidity,  temperature  and  pressure,  (U,  t 
and  P).    Thus  one  way  of  assuring  that  the  total  uncertainty  of      is  within  accept- 
able limits  is  to  make  sure  that  the  individual  uncertainties  from  U,  t  and  P  are 
each  less  than  1//3  of  the  total  acceptable  uncertainty.    That  is,  we  set  a  goal  that  . 
the  contribution  to  uncertainty  in  p^  from  U,  P,  or  t  should  be  less  than  0.0017  mg/cm' 

We  will  see  in  Section  2.4  that  this  goal  can  be  met  or  exceeded  with  a  modest  invest- 
ment  in  equipment.    We  now  show  the  relative  importance  of  measurements  of  humidity, 

3 

The  CCM  recommends  that  this  equation  referred  to  as  "Euqation  for  the  determination 
of  the  density  of  moist  air  (1981):  in  anticipation  of  improved  reference  data.  There 
are  some  stylistic  differences  in  the  CCM  version  of  [2]  but  the  calculated  values  are 
completely  consistent  at  meaningful  levels  of  precision. 

Equations  (20a)  and  (20b)  are  somewhat  simplified  (see  Appendix  B).    Weighing  errors 
due  to  the  simplifications  are  usually  very  small  and  are  estimated  in  Table  5,  p.  427 
of  ref.  [2]. 
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temperature,  and  pressure  to  the  uncertainty  in  air  density;  we  also  indicate  now 
accurately  these  must  be  measured  to  achieve  our  goal  of  0.0017  mg/cm3  for  each  one. 

2.3.1.  Humidity 

It  is  easiest  to  discuss  the  effects  of  humidity  (U)  via  examples.  For  practi- 
cal purposes  we  can  consider  as  standard  a  condition  when  t  =  20  °C,  P  =  760.00  mm  H< 
and  U  =  50  percent. 

Example  1 

Consider  t  =  20  °C,  P  =  760.00  mm  Hg,  U  =  50  percent  (i.e.  standard  con- 
dition).   Using  eq.  (20b),  we  find 

PA  =  1.2045  -  0.0053  =  1.1992  mg/cm3  . 

3 

We  note  that  the  humidity  correction,  0.0053  mg/cm  ,  corresponds  to  a 
correction  of  only  0.44  percent  in  p^. 

Example  2 

Same  conditions  as  in  Example  1  except  t  =  30  °C.    We  calculate 

PA  =  1.1648  -  0.0093  =  1.1555  mg/cm3 
and  the  humidity  corrections  is  ^0.8  percent  of  p^. 
Example  3 

Same  as  Example  1  except  U  =  60  percent. 

PA  =  1.2045  -  0.006330  =  1.1982  mg/cm3 

so  that  the  humidity  correction  is  0.5  percent. 

Summarizing  the  results  on  humidity  we  can  make  the  following  observations: 

i.     The  magnitude  of  the  humidity  correction  to  the  air  density  at  standard 

conditions  is  ^0.4  percent. 

ii.     By  increasing  the  temperature  from  20  °C  to  30  °C,  the  humidity  correction 

amounts  to  0.8  percent  of  the  air  density.    Thus  even  under  these  high 

temperatures,  very  unlikely  in  a  temperature-controlled  laboratory,  the 

humidity  correction  term  is  still  quite  small. 

iii.     A  change  of  10  percent  in  the  relative  humidity  (i.e.  U  =  Un+10)  changes 

•    J           u       i    n  no           +    t>       1.1992-1.1982N  u 
the  air  density  by  only  0.08  percent    (i.e.  — ^  -j-g-gg  ). 

We  can  conclude,  therefore,  that  if  we  want  to  know  the  air  density  to  within 

0.0017  mg/cm    then  the  relative  humidity  should  be  known  to  ±16  percent.    (If  a 
safety  factor  of  4  is  required  then  the  relative  humidity  should  be  known  at  the  ±4 
percent  level ) . 
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2.3.2.  Temperature 

The  effects  of  temperature  on  the  humidity  term  were  already  demonstrated  in 
Example  2.    In  particular  we  note  that  a  10  °C  change,  yielded  a  humidity  correction 
of  0.8  percent. 

The  most  sensitive  effect  of  temperature  on  calculated  air  density,  however, 
occurs  in  the  first  term 

0.46460  D 
PA  =    P 

rt  (t+273.15) 

From  Example  1,  we  note  that  at  standard  conditions  the  first  term  yields 

p  =  1.2045  mg/cm3 

Example  4 

P  =  760.0  mrn  Hg,  t  =  21  °C 

In  this  case  the  first  term  of  eq.  (20b)  gives 

.3 


PA  =  1.2004  mg/ 


cm 


3 

and  pA  has  changed  from  p  by  0.0041  mg/cm   or  0.34  percent. 

3 

To  achieve  an  uncertainty  of  0.0017  mg/cm  ,  we  must  measure  the  temperature 
accurately  to  ^±0.4  °C  (which  reduces  to  ±0.1  °C  for  a  safety  factor  of  4.) 

2.3.3.  Pressure 

The  arguments  for  the  pressure  uncertainty  follow  those  of  the  temperature. 
Once  again  we  only  have  to  consider  the  first  term  of  the  density  equation 

0.46460  q 
p    =  p 

M  (t+273.15) 

For  standard  conditions  of  t  =  20  °C,  and  P  =  760.0  mm  Hg 


P  =  1.2045  mg/cm3 


Example  5 


P  =  764.00  mm  Hg 
t  =  20  °C 

In  this  case  pA  =  1.2108  and  has  changed  from  p  by  0.53  percent. 

Thus  in  order  to  know  pA  to  0.0017  mg/cm3  we  will  need  to  have  accurate  baro- 
metric measurements  to  ±1.1  mm  Hg  (which  reduces  to  ±0.3  mm  Hg  for  a  safety  factor  of 

2.3.4.    Overall  Uncertainty  of  Air  Density 

We  have  seen  in  eq.  (20)  that  the  air  density  is  a  function  of  the  temperature, 

pressure,  and  humidity 
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To  evaluate  the  overall  uncertainty  of  p«,  we  assume  uncorrelated  errors  and  develop 
a  final  value  in  Appendix  C: 

6pA 

 0.25  percent  . 

PA 

Note  in  particular  that  this  value  was  obtained  using  the  errors  established  in 
Sections  2.3.1,  2.3.2,  and  2.3.3  of  this  section  (with  no  safety  factor). 

2.3.5.    Effects  of  Small  Changes  in  Air  Density  on  the  Mass  Value 

Having  established  the  sensitivities  of  U,  t,  P  on  p^,  the  air  density,  we  can 
now  proceed  and  see  what  effect  a  small  change  of  has  on  the  mass  determination. 
The  arguments  can  best  be  carried  through  with  the  following  example: 

Assume  that  we  are  comparing  the  masses  of  two  1  kg  weights. 

3 

Weight  #1  has  density  of  8.4  g/cm 

3 

Weight  #2  has  density  of  7.5  g/cm 

3 

Air  density  p^  =  1.20  mg/cm 

If  the  air  density  is  known  to  an  accuracy  of  0.25  percent,  what  uncertainty  does 
this  cause  in  the  mass  difference  calculations? 

Answer:     We  found  in  Section  2.1.5  that  the  uncertainty  for  true  as  well  as  apparent 
mass  was  given  as 

uncertainty  =  ±<5pA  (Vqi~V02^ 

=  ±6pAN  (1-  -  ]-) 
pl  p2 

where  we  have  approximated  the  volumes  by 

V]  ~  N/pp  V2  *  N/p2 

with  N  being  the  nominal  mass  of  each  weight. 

Substituting  numerical  values,  yields 

Uncertainty  =  ±0.0025  x  p«  x  1000  (—  — ) 

A  8.4  7.5 

=  ±0.0025  x  1.2  x  1000  (-  1 


8.4  7.5 
=  ±0.043  mg. 
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This  is  about  one-half  the  calibration  uncertainty  limit  provided  by  the  NBS  for 

3 

high-quality  1 -kg  weights.    The  span  7.5  to  8.4  g/cm    includes  the  densities  of  1 -kg 
weights  used  in  routine  fine  work. 

It  is  usually  true  that  the  most  stringent  requirements  placed  on  buoyancy 
measurements  are  in  the  precise  comparison  of  kilograms.    An  uncertainty  in  p«  of 
0.25  percent  is  usually  sufficient  even  in  this  case. 

CAUTION:    Some  weights  are  manufactured  for  purposes  other  than  general  mass 
metrology.    For  instance,  a  scientist  may  have  a  1 -kg  aluminum  weight  which  needs  the 
best  calibration  possible.    In  such  special  situations,  the  requirements  and  limita- 
tions placed  on  buoyancy  work  must  be  examined  on  a  case-by-case  basis. 

2.4.    Measuring  Equipment  for  Air  Density 

The  previous  section  yielded  information  on  the  accuracies  required  for  the 
three  parameters  needed  for  the  air  density.    In  this  section  we  discuss  briefly  the 
types  of  equipment  required  to  meet  our  accuracy  goals. 

2.4.1.  Humidity 

This  parameter  is  usually  measured  with  a  psychrometer  or  hygrometer.  Several 
commercial  instruments  can  provide  accuracies  of  ±2%  relative  humidity  over  a  wide 
range  such  as  10-80%.    The  calibration  of  such  instruments  can  be  verified  at  fixed 
points  by  the  user  by  means  of  an  in-laboratory  calibration  station.    Since  an 
accuracy  of  ±16%  (or  ±4%  using  a  "safety  factor  of  4")  is  required  for  mass  metrology, 
these  instruments  are  more  than  adequate. 

2.4.2.  Temperature 

An  accuracy  of  ±0.4  °C  (or  ±0.1  °C  for  a  "safety  factor  of  4")  is  required. 
Simple,  mercury-filled,  glass  thermometers  can  provide  reading  accuracies  up  to  ±0.1 
°C.    More  elaborate  units  can  yield  ±0.01  °C. 

The  user  should  be  very  careful  in  taking  temperature  readings.    Generally,  the 
physical  presence  of  the  metrologist  tends  to  warm  up  the  balance  with  respect  to  the 
rest  of  the  room.    It  is  therefore  essential  that,  no  matter  what  form  of  thermometer 
is  used,  its  sensor  be  placed  as  close  as  possible  to  the  balance  pan.    Also,  every 
attempt  should  be  made  to  ensure  that  objects  being  weighed  are  in  thermal  equilib- 
rium with  the  balance.    The  latter  requirement  is  especially  important  when  comparing 
weights  of  large  surface  area  or  of  different  geometries. 

2.4.3.  Barometric  Pressure 

The  required  accuracy  is  ±1.1  mm  Hg  (or  ±0.3  mm  Hg  for  a  "safety  factor  of  4"). 
One  can  obtain  mercury  manometers  or  aneroid  barometers  with  reading  accuracies  of 
±0.05  mm  Hg.    Since  this  is  sufficient  for  usual  mass  calibrations,  one  does  not  have 
to  push  for  higher  accuracy  instruments. 

Present  technology  can  provide  defining  instruments  with  overall  uncertainty  of 
±0.01  mm  Hg.    Since  such  devices  are  the  most  accurate  on  the  market,  they  are 
generally  used  in  calibration  of  other  manometers  or  aneroid  barometers  which  con- 
sequently would  have  greater  uncertainty. 
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Many  laboratories  possess  either  Fortin-type  or  aneroid  barometers.  These 
require  calibration  against  some  defining  instrument—such  as  a  mercury  manometer  or 
a  piston  gage;  but  such  defining  instruments  are  costly  and  may  not  be  available  in 
the  laboratory.    On  the  other  hand,  aneroid  and  Fortin-type  barometers  are  best 
calibrated  in  place. 

A  solution  to  this  problem  has  recently  been  demonstrated  [9].    Two  weights  of 
well-known  and  nearly-equal  mass  but  having  very  different  volumes  are  compared  in 
the  laboratory  on  a  sensitive  balance.    The  measured  difference  in  balance  readings 
between  these  two  weights,  kA6,  determines  the  air  density  in  the  balance  enclosure 
via  eq.  (9).    The  air  density  is  also  calculated  from  eq.  (20)  using  a  calibrated 
thermometer,  a  calibrated  hygrometer,  and  the  uncalibrated  barometer.    The  difference 
in  the  two  measurements  of  air  density  serves  to  calibrate  the  barometer.  Calibration 
uncertainties  of  less  than  1  mm  Hg  (three  standard  deviations  plus  known  systematic 
errors)  have  been  demonstrated. 

In  summary,  we  note  that  all  the  requirements  for  air  density  measurements  can 
be  met  with  presently  available  instrumentation.    If,  however,  the  accuracy  of  ±0.25 
percent  for      has  to  be  improved,  then  the  limiting  instrument  will  be  the  pressure 
gage.    (It  does  not  make  sense  to  improve  temperature  and  humidity  resolution  by  an 
order  of  magnitude,  if  the  pressure  cannot  be  read  with  higher  precision.) 
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3.    WEIGHING  METHODS 


Until  this  point,  we  have  assumed  that  all  mass  comparisons  have  been  carried 
out  on  the  idealized  balance  defined  in  Section  2.    In  this  section  we  develop  the 
means  of  evaluating  the  constant,  k  of  eq.  (8).    We  also  extend  our  analysis  to 
actual  balances  currently  in  use.    Just  as  one  need  not  be  a  mechanical  engineer  to 
drive  a  car,  so  one  need  not  have  a  detailed  knowledge  of  balance  design  in  order  to 
carry  out  successful  mass  measurements.    Nevertheless,  a  general  knowledge  of  balance 
design  is  often  useful.    The  review  paper  by  Schoonover  [10]  provides  an  excellent 
introduction  to  this  subject. 

The  idealized  balance  of  Section  2  had  no  imprecision.    Unfortunately,  this  is 
not  the  case  for  real  balances.    By  saying  a  balance  has  imprecision,  we  mean  two 
distinct  contributions  to  uncertainty: 

1)  Resolution.    A  balance  observation  can  be  no  more  certain  than  one's 
ability  to  resolve  the  least-significant  digit  of  the  read-out. 

2)  Reproducibility.    A  balance  observation  can  be  no  more  certain  than  the 
tendency  of  the  balance  to  produce  identical  read-outs  under  identical 
conditions.    The  lack  of  this  property  is  called  "scatter." 

For  balances  used  in  the  best  metrology,  scatter  is  by  far  the  more  important  of 
the  two.    What  we  loosely  term  the  scatter  of  the  balance  is  not  entirely  intrinsic 
to  the  balance.    That  is,  it  may  depend  on  the  type  of  table  on  which  the  balance  is 
placed,  the  type  of  air-conditioning  in  the  laboratory,  the  skill  of  the  balance 
operator,  etc.    We  need  a  measure  of  the  scatter  which  includes  all  these  effects. 
Such  a  measure  is  provided  by  the  "process  standard  deviation"  [7,8]. 

The  standard  deviation  for  a  particular  weighing  process  is  easily  estimated  by 
repeating  the  process  a  number  of  times.    (Ten  repetitions  are  usually  sufficient  to 
give  a  respectable  estimate.)    Two  weights,  representative  of  the  weights  normally 
measured  on  a  particular  balance,  are  intercompared  n  times  by  one  of  the  methods 
described  below.    The  n  values  obtained  are  aq ,  a9,  aq,         a    -.,  a  .    The  average 
of  the  n  values  is  i     ^     j  n-i  n 


n 

.S  a. 
a  =  i=L_L 

n 

and  the  standard  deviation  (a)  of  the  balance  is  estimated  to  be 

n  \  1/2 

s  =   3   I, 

n-1 

with  s  ■+  a  as  n  «>.  (The  estimate  s  approaches  the  actual  standard  deviation  a  as 
n  is  made  larger.) 

The  standard  deviation  should  be  estimated  in  this  way  at  regular  intervals  in 
order  to  spot  trends  or  shifts  in  the  evolution  of  s  with  time.    A  slow  rise  in  s 
could  signal,  for  example,  a  slow  deterioration  of  the  balance.    An  abrupt  shift 
could  signal  a  problem  with  the  balance,  a  change  in  the  air-conditioning  system  of 
the  laboratory,  or  some  other  problem  which  requires  attention. 
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Before  actually  discussing  various  types  of  weighing  procedures,  we  have  to 
define  some  terms  required  in  the  descriptions.    The  following  symbols  will  be 
util ized. 

M  =   any  unknown  weight  with  true  mass  MT  whose  value  has  to  be  determined 

C  =   a  counterpoise  (or  counterweight)  with  true  mass  CT. 

This  weight  is  usually  required  on  beam  balances.    It  can  be  of  any  stable 
shape,  form  and  composition.    However,  it  has  to  be  adjusted  so  that  the 
scale  (balance)  does  not  deflect  past  its  maximum  allowable  value  when 
loaded  with  a  weight.    This  mass  appears  in  equations  which  are  combined  in 

such  a  way  that  CT  drops  out  of  the  final  result. 
S  =   a  standard  weight  with  true  mass  S  . 

This  is  a  known  mass  against  which  the  unknown  is  measured.    It  is  always 

assumed  that  in  any  process  ST  is  known  and  is  traceable  to  the  NBS.  In 
the  following  description  S  can  be  of  a  single  unit  or  can  be  made  up  of 


S    =  Z    S.  such  that  the  sum  S    roughly  equals  M   or  C  . 

i=l  1 

A  =   a  sensitivity  weight  with  true  mass 

This  mass  is  usually  of  very  small  magnitude.    In  principle  it  is  only  used 
to  measure  the  deflection  of  the  scale  per  unit  mass.    This  mass  is  very 
critical.    It  must  be  well  calibrated  and  traceable  to  the  NBS,  although 
such  calibration  is  a  relatively  simple  matter. 

Of  the  various  weighing  procedures,  we  will  mention  Direct  Weighing,  Direct- 
Reading  Weighing,  and  discuss  Substitution  and  Transposition  Methods. 

3.1 .    Direct  Weighing 

This  type  of  weighing  is  seldom  used  in  the  mass  calibration  program.  It 
requires  the  use  of  an  equal-arm  balance.    Essentially,  the  reading  steps  are  as 
follows: 

1.  Release  brake  and  take  null  reading  with  both  pans  empty. 

2.  Load  one  pan  with  M  and  balance  with  S  masses  on  the  opposing  pan  until  the 
null  position  is  once  again  reached. 

We  have  as  final  result,  if  buoyancy  corrections  are  insignificant, 

T      n  T 

M    =    E    S. 1 

i=l  1 
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3.2.  Direct-Reading  Weighing 

This  weighing  procedure  requires  a  one-pan  balance  that  reads  directly  in  S.I. 
mass  units  or  some  other  commercially  recognized  unit  (e.g.  pounds  or  carats).  The 
procedural  steps  are  as  follows: 

1.  Null  the  balance  without  any  load  on  the  pan. 

2.  Place  M  on  pan  and  record  the  mass  directly. 

This  method  is  obviously  quite  simple  but  usually  is  not  accurate  enough  for  metrol- 
ogy.   It  is  appropriate  for  scales  most  commonly  used  in  stores.    Nevertheless,  the 
method  is  quite  similar  in  principle  to  that  described  in  Section  3.4.5.,  below. 

3.3.  Substitution  Weighing 

This  weighing  procedure  is  the  usual  one  employed  for  both  single  pan  or  equal  - 
arm  balances.    The  basic  idea  is  the  comparison  of  an  unknown  mass  with  a  standard 
mass  and  a  sensitivity  mass. 

3.3.1.    Single  Pan  Balance 

Known  as  "direct-reading  analytical  balances,"  they  contain  built-in  weights 
which  are  usually  manipulated  by  external  knobs  or  dials.    These  we  will  refer  to  as 
"dial  weights."    The  remaining  balance  reading  is  indicated  by  the  rest  point  of  a 
moving  optical  screen  or  by  an  automatic  electronic  display.    This  part  of  the 
balance  output  is  called  the  "screen  reading."    Thus  each  balance  indication  is  the 
summation  of  the  dial  weights  and  the  screen  reading.    Calibration  of  the  dial 
weights  is  unnecessary  for  the  measurements  described  below.    The  majority  of  balance 
operations  outside  metrology  (e.g.,  chemistry,  metallurgy,  etc.)  do  rely  on  the 
accuracy  of  the  built-in  dial  weights,  however.    For  calibration  of  these  weights  the 
reader  is  referred  to  Ref.  [11]  and  Appendix  A. 

3.3.1.1.    Single  Substitution 

To  find  the  mass  difference  between  two  weights  of  nominally  equal  value,  we 
will  first  describe  the  method  of  "single  substitution".    For  the  moment,  let  us 
assume  that  the  weights  are  matched  closely  enough  so  that,  when  placed  on  the 
balance,  they  both  require  the  same  setting  of  dial  weights.    Therefore  they  differ 
only  in  their  screen  readings.    If  we  look  into  the  balance  (fig.  3),  we  see  that  the 
pan  is  suspended  from  a  beam  which  pivots  about  a  very  good  bearing  (the  fulcrum, 
usually  a  knife  on  a  flat).    On  the  end  of  the  beam  opposite  to  the  pan,  is  a  fixed 
counterpoise  C.    (There  will  also  be  either  a  servo-motor  on  this  end  or  a  dashpot  to 
damp  the  beam  oscillations.    These  details  are  unimportant  to  our  derivation.)  Also 
suspended  from  the  pan-end  of  the  beam  are  the  built-in  dial  weights.    Thus  when  an 
object  is  placed  on  the  pan,  dial  weights  are  removed  so  that  the  force  of  the  total 
assembly  on  the  pan  side  (object  being  weighed  plus  remaining  dial  weights)  balances 
the  counterweight  as  closely  as  possible.    The  remaining  imbalance  causes  the  beam  to 
rock  slightly  out  of  its  horizontal  position.    The  angle  of  tilt  is  proportional  to 
the  remaining  force  imbalance.    In  a  mechanical  balance,    this  angle  is  directly 
proportional  to  the  screen  readings.    In  a  servo-controlled  balance,  a  motor  is  used 
to  drive  the  beam  back  to  horizontal.    The  screen  readings  are  then  proportional  to 
the  average  force  generated  by  the  motor. 
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In  the  balances  just  described,  we  can  define  the  following  quantities: 
ZdT  =  the  true  mass  of  the  dial  weights  removed 
EV.J  =  the  volume  of  the  dial  weights  removed 
6-j    =  the  screen  reading 

6q    =  the  screen  reading  when  the  pan  is  empty  and  the  dial  weights  are  replaced 
(this  was  always  zero  for  the  ideal  balance  of  Section  2) 

k     =  the  proportionality  constant  between  screen  reading  and  mass. 

(Recall  the  previous  definitions  for  M,  S,  A.)    We  also  recall  from  eq.  (8)  that 
for  two  weights,       and  1^  with  volumes  V-|  and  V^,  the  weighing  equation  can  be 

written  as  (using  true  masses) 

g  (M^-PrnV  =  9  (M2T-pmV2}  +  9  k<Ae)  ■ 

With 

Pm  =  PA  =  air  density; 
(A6)  =  e1  -  e2; 

and  setting 

=  edT 
v2  =  zv, 

we  arrive  at^ 

g(M{  -  paV^  =  (edT  -  pA  EV.)g  +  gk(61  -  eQ) 

where  6q  e      in  this  notation. 

The  following  observations  can  be  made: 

t   The  counterpoise  weight  does  not_enter  into  the  final  result. 

•   The  gravitational  acceleration,  g,  can  be  cancelled  from  both  sides  of  the 

above  balance  equation, 
t   The  above  equation  is  identical  to  that  for  the  true  mass  difference 

between  weight  #1  and  the  removed  dial  weights  (see  Appendix  A). 

In  single  substitution,  we  compare  a  weight  M  with  a  standard  weight  S.    Either  M 
or  S  may  actually  be  a  summation  of  weight  pieces.    A  small  calibration  weight, 
sensitivity  weight  A,  is  also  required. 


A  full  derivation  would  involve  a  balance  equation  for  each  weighing  operation  (i.e. 
the  pan  unloaded  and  the  pan  loaded).    These  equations  would  each  involve  the  counter- 
poise C.    Combining  these  relations  would  then  give  the  desired  result. 
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The  steps  in  single  substitution  are: 

1)  Load  M  and  record  screen  reading  6-. 

2)  Remove  M 

3)  Load  S  and  record  screen  reading  6£ 

4)  Add  A  to  S  and  record  screen  reading  63 

In  general  we  adjust  S  and  choose  A  so  that 

•  The  same  dial  weights  are  used  in  steps  1,  3, and  4. 
t  |e3  -  62|is  at  least  four  times  as  large  as  |e2  - 

•  (63  -  62)  ^  ^  *he  total  range  of  the  screen. 

From  the  above  steps  1,  3  and  4,  we  have  the  following  relations: 

ke1  =  (mT-paVm)  -  (zd{-paev1 )  +  ke0  (22) 
ke2  =  (st-pavs)  -  (zdT-paev.)  +  keQ  t  (23) 

ke3  =  (st-pavs)  +  (at-pava)  -  (zdT-paev.)  +  ke0  •  (24) 

The  symbol  =  indicates  that  we  are  now  assuming  that  some  or  all  of  the  measured 
values  have  a  non-negligible  uncertainty. 

Subtracting  (23)  from  (22)  yields 

mt  -  sT  =  pa(vm-vs)  +  k(ere2)    .  (25) 

Subtracting  (24)  from  (23)  yields 

k  =  7  H   •  (26) 


(e3"e2)  ' 


Finally,  substituting  for  k  from  (26)  into  (25)  yields 

mt-st  =  pa(vm-vs)  +    i  Pk\  (ere2)   .  (27) 


(e3-62) 


Note  that  this  equation  is  identical  to  (10),  the  crucial  difference  being  that 
we  have  found  the  value  for  the  proportionality  constant  k. 

One  may  wonder  why  k  must  be  evaluated  for  every  measurement  if  it  is  a  con- 
stant.   The  answer  is  that  for  most  sensitive  balances,  k  is  truly  constant  only  for 
short  periods  of  time  and  only  over  reduced  regions  of  the  range  of  the  screen. 

In  the  derivation  above  we  choose  A  so  that  63-62,  is  about  20  percent  of  the 

screen  range  (fig.  4).    In  this  reduced  region,  k  will  be  sufficiently  constant  for 
metrology  if  one  is  using  well-made  balances. 
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FORCE  ON  BALANCE 


Figure  4.    Balance  indications  as  a  function  of  the  force  on  the  pan.    The  response 
shown  is  that  of  the  screen  when  the  dial  weights  are  not  changed. 
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3.3.1.2  Double-Substitution 

In  the  derivation  of  single-pan  substitution  given  above,  we  have  treated  6Q  as 
a  constant.    This  is  the  same  as  assuming  that  the  balance  indication  does  not  drift 
with  time.    But  this  is  rarely  the  case.    Usually,  however,  the  zero-drift  in  the 
balance  will  be  approximately  linear  over  a  short  period  of  time.    The  method  known 
as  single-pan  double-substituion  has  the  following  two  advantages  over  single- 
substitution: 

•  linear  drifts  in  the  balance  zero  are  eliminated  if  there  is  an  equal  time 
between  successive  observations; 

•  a  second  estimate  of  the  mass  difference  between  M  and  S  is  obtained. 

Both  benefits  are  obtained  by  a  single  additional  measurement: 

5)  Remove  S  from  pan,  leave  A  where  it  is 

6)  Place  M  on  the  pan  and  record  secreen  reading  6*.    The  last  step  gives  the 
additional  balance  relation 

ke4  =  (mt-pavm)  +  (at-pava)  -  (ed]"-paev.)  +  ke0   .  (28) 

Substracting  (24)  from  (28)  yields 

mt  -  sT  =  pa(vm-vs)  +  k(e4-e3)  . 

Substituting  for  k  from  eq  (26)  then  gives 

T      T  (aT-pAVA} 

PA<  W  +  -77-77-    (Ve3>    '  <29> 
^3  2' 

Finally,  by  adding  (29)  and  (27)  we  obtain 

pa(vh-vs)  ♦  — (ere2-e3+94)   .  (30) 

*  3"  2' 

One  might  wonder  how  eq.  (30)  would  have  to  be  changed  if  6Q  were  drifting 

linearly  with  time  instead  of  being  constant.    The  answer  is:    the  term  (0-| "e2_e3+e4^ 

would  still  be  strictly  correct;  the  term  (63-62)  would  not  be  strictly  correct. 

However,  provided  that  the  drift  among  readings  is  small  compared  to  (63-62),  the 

error  from  this  source  will  be  negligible.    This  will  in  actuality  always  be  the 
case.    If  it  were  not,  the  balance  would  be  drifting  so  badly  that  no  stable  readings 
could  be  obtained.    If  instead  of  (63-62)  one  substitutes  the  term  (363-3e2+e-j -e4)/2 , 

linear  drift  in  the  sensitivity  calculation  will  be  removed  [12].    This  more  com- 
plicated formulation  is  generally  used  when  data  are  being  analyzed  by  digital  - 
computer  software. 
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CAUTION:    If  6j  or  62  is  very  close  to  zero,  balance  drift  may  cause  readings  to 

become  negative  or  to  go  off  scale.    To  avoid  this  problem,  a  small,  uncalibrated 
"tare  weight"  can  be  added  to  the  pan  during  all  measurements.    Mathematically,  this 
is  equivalent  to  changing  6q  by  a  constant  amount  and  so  has  no  effect  on  the  final 
results. 

3.3.2.    Two-Pan  Equal -Arm  Balance 

This  is  the  oldest  type  of  precision  balance.    Its  basic  design  is  represented 
by  the  "scales  of  justice"  or  by  the  zodiac  sign  for  Libra. 

3.3.2.1.  Single  Substitution 

This  method  is  very  similar  to  single-substitution  on  a  single-pan  balance.  The 
read-out  provided  by  an  equal -arm  balance  corresponds  to  just  the  screen-reading  of 
the  single-pan  balance.    An  uncalibrated  counterpoise  which  has  nearly  the  same  mass 
as  M  and  S  is  placed  on  one  of  the  pans  throughout  the  measurement  sequence.  The 
size  of  the  counterpoise  is  selected  to  ensure  that  the  balance  indication  will  be  on 
scale  with  M  on  the  second  pan.    The  weights  M,  S  and  a  sensitivity  weight  A  are 
placed  on  the  second  pan  in  nearly  the  same  sequence  used  in  the  one-pan  case. 

The  required  steps  are  as  follows: 

1)  Place  C  on  first  pan,  M  on  second  pan  and  record  indicator  reading  6-j . 

2)  Remove  M 

3)  Place  S  on  second  pan  and  record  indicator  reading         (Recall  that  S  is 
a  summation  of  calibrated  weights  which  is  chosen  to  be  near  the  scale 
indication  of  M. ) 

4)  Place  the  sensitivity  weight  A  on  that  pan  which  will  cause  the  indicator 
to  move  toward  the  center  of  the  scale.    Record  indicator  reading  6g. 

It  is  convenient  to  choose  as  the  second  pan  (i.e.  the  one  on  which  M  and  S  are 
placed)  the  one  for  which  increasing  increments  of  load  give  increasingly  positive 
indicator  readings. 

Going  through  detailed  arguments  similar  to  those  for  single-pan,  single  substi- 
tution, we  arrive  at 

mT-  s  =         +  ir-rr  (ere2>  (31) 

|e3"e2| 

The  absolute  value  j^-^l  comes  about  because  the  A  weight  might  be  placed  on 

either  of  the  two  pans.    Note  that  if  the  arms  of  the  balance  are  not  equal,  the 
sensitivity  weight  must  be  placed  on  the  second  pan. 

3.3.2.2.  Double-Substitution 

As  in  the  case  of  the  single-pan  balance,  one  additional  reading  is  taken  for 
double  substitution. 

5)  Remove  S  from  the  balance;  leave  A  where  it  is, 

6)  Place  M  on  the  second  pan  and  record  the  indicator  reading  8>. 
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By  obvious  extension  of  the  preceding  arguments,  we  arrive  at  the  equation  for 
double-substitution  on  an  equal -arm  balance: 

MT  -  ST  =  PA(VM-VS)  +  -r— M       1    2    3   4       .  (32) 

1 0  2 —®  2 1 

Again,  to  be  strictly  correct  in  the  presence  of  linear  drfit,  the  term   63-62  can 
be  replaced  by    (3e3-362+e^-e4)/2.    Double  substitution  in  this  case  has  the  same 
advantages  as  weighing  on  a  single-pan  balance. 

In  summary,  we  have  derived  the  following  equations  for  the  four  types  of 
substitution  weighing  we  have  discussed: 

Balance  Type  Substitution  Type  Equation 

Single-Pan  Single  (27) 

Single-Pan  Double  (30) 

Two-Pan  Single  (31) 

Two-Pan  Double  (32) 

3.4.    Transposition  Weighing 

Two-pan,  equal -arm  balances  can  also  be  used  to  do  transposition  weighing. 
This  method  does  not  require  a  separate  counterpoise  and  also  doubles  the  balance 
sensitivity. 

3.4.1.    Single  Transposition 

The  following  steps  are  required: 

1.  Place  N!  on  one  of  the  pans,  e.g.,  the  left  pan  and  S  on  the  right  pan.  S 
should  be  adjusted  so  that  the  indicator  reads  on  scale.    Record  the  indication 

er 

2.  Remove  M  and  S. 

3.  Replace  M  and  S  on  the  balance  but  in  transposed  position  (that  is,  M  on  the 
right  pan  and  S  on  the  left  pan).    Record  the  reading 

4.  Add  A  to  that  pan  which  will  cause  the  indicator  to  deflect  towards  the  center 
of  the  reading  scale.    Record  the  indication  63. 

The  three  indications  give  the  relations: 


MT  -  sT  =  pa(vm-vs)  +  k(ere0)  (33) 

sT  -  mt  =  pa(vs-vm)  +  k(e2-e0)  (34) 

and 

sT  -  mt  =  pa(vs-vm)  ±  (at+pava)  +  k(e3-e0)  .  (35) 

Here  again,  6q  is  the  balance  indication  when  both  pans  are  empty.    The  ±  sign 

is  used  to  indicate  that,  in  the  last  equation,  the  sign  of  (AT-pAVA)  depends  on 
which  pan  held  the  sensitivity  weight. 
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Subtracting  (34)  from  (35)  yields 

J. 

k  -  t  =-£    •  (36) 

'3-92[ 

where,  once  again,  the  absolute  value  was  used  to  take  account  of  whichever  pan  the 
sensitivity  weight  is  placed  on. 

By  subtracting  (34)  from  (33),  we  have 

mt  -  sT  =  pa%vs)  +  |  (ere2)   .  (37) 

Finally,  we  can  substitute  from  (36)  into  (37)  to  obtain 


MT    J  -    rv  v  \  ,  (aT'paVa}  (ere2} 

pa(vm-vs)  +  ~—   . 


'3"D2 


(38) 


CAUTION:    In  transposition  weighing,  attention  has  to  be  paid  to  the  sense  of 
the  balance  scale.    We  have  assumed  above  that  the  balance  indication  becomes  more 
positive  as  the  mass  on  the  left  pan  is  increased.    For  scales  of  opposite  sense, 
the  plus  sign  of  eq.(38)  should  be  changed  to  minus.    Also,  if  the  arms  of  the 
balance  are  insufficiently  "equal,"  transposition  weighing  will  result  in  error. 

3.4.2.    Double  Transposition 

The  procedure  for  this  method  follows  that  outlined  above  but  with  the  follow- 
ing additions: 

5)  Remove  M  and  S,  leaving  A  where  it  is. 

6)  Place  M  on  the  left  pan  and  S  on  the  right  pan.    Record  indication  0*. 

These  steps  provide  us  with  the  additional  equation 

mt  -  sT  =  pa(vm-vs)  ±  (at-pava)  +  k(e4-e0)  (39) 

where  the  ±  sign  is  needed  once  again  since  the  location  of  the  A  weight  is  ambigu- 
ous . 

By  subtracting  (35)  from  (39)  we  have 

mt  -  sT  =  pA(vM-vs)  +  |  (e4-e3) 

and  substituting  for  k  from  (36)  yields 


T       T  1  (AT-pAVA} 

d  e3-e2 
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Adding  (40)  to  (38)  finally  results  in 


t      t  IaT-PaVa 

m  -  s  =  pa(vm-vs)  +  \ — (ere2-e3+e4)   .  (41) 


4 


3  u2 


Again  the  reader  is  cautioned  that:  the  plus  sign  in  front  of  the  second  term  applies 
only  to  balances  whose  indications  become  more  positive  as  the  mass  in  the  left  pan 
is  increased;  in  the  presence  of  linear  drift,  the  term   6^-99   can  be  replaced  by 
(363-362+0^)72.  J  c 

3.5.    Weighing  on  Electronic  Balances 

The  most  recent  development  in  mass  technology  is  the  appearance  of  the  com- 
pletely electronic  balance.    This  can  be  thought  of  as  a  single-pan  balance  whose 
screen  span  has  become  so  large  that  dial  weights  are  no  longer  necessary.  These 
balances  are  rugged,  easy  to  use,  and  most  are  directly  interfaceable  to  computers  or 
data  loggers.    Unfortunately  at  present  they  rarely  attain  the  precision  achieved  by 
the  best  mechanical  balances.    For  metrology  purposes,  these  balances  are  used  just 
as  the  single-pan  balances  discussed  above  and  the  same  equations  apply.    If  and  only 
if  experience  verifies  that  the  constant  k  is  indeed  constant,  can  weighing  pro- 
cedures be  simplified  so  that  the  sensitivity  weight  need  not  be  measured  for  every 
substitution  or  double  substitution. 

The  results  we  have  just  derived  can  be  generalized.    First,  although  only  true 
mass  has  been  dealt  with,  the  reader  should  by  now  be  able  to  generate  the  corre- 
sponding equations  for  apparent  mass.    Next,  since  S  is  assumed  to  be  an  assembly  of 
calibrated  weights,  all  the  measurements  described  in  this  section  provide  an  esti- 
mate for  M  ,  although  they  do  not  constitute  an  actual  calibration. 

Recall  that  we  had  (eq.  (10)) 

MT  =  ST  +  pA  (VM-VS)  +  kA6  (10') 

or  equivalently 

T  ST-pflVQ+kA6 
MT  =   £-S   .  (10") 

1 -PA/p 

Usually,  the  manufacturer  of  a  weight  supplies  the  buyer  with  the  weight's  density  so 
that  formulation  (10")  is  more  appropriate.  Equation  (10 l)  is  often  used  in  computer 
analysis  of  data,  however.    In  this  case,  one  estimates  VM  as  VM     N/p,  where  N  is 

the  nominal  mass  value  of  M.    Equation  (10')  then  gives  a  first  order  approximation 

for  MT.    This  first  estimate  can  be  used  to  find  a  better  value  for  VM  and  the 
process  repeated  until  convergence  is  reached. 

If  S  and  M  are  both  unknown  and  not  sufficiently  close  in  mass,  known  weights 
can  be  added  to  either  S  or  M.    These  added  weights  are  then  included  in  the  weighing 
equations. 
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Example 

Sometimes  M  itself  will  be  a  summation  of  weights.    Let  us  say,  to  take  a 
specific  example,  that  M  is  the  sum  of  three  weights  whose  nominal  values  are  50  mg, 
30  mg,  and  20  mg.    Furthermore,  the  weight  manufacturer  specifies  that  the  50-mg  and 

3  6 

30-mg  weights  are  made  of  tantalum  (p  =  16.6  g/cm  ;  a  =  20  x  10    /°C)  but  the  20-mg 

3  6 

weight  is  made  of  aluminum  (p  =  2.7  g/cm  ;  a  =  69  x  10"  /°C).  What  is  the  effective 
density  of  the  summation  and  what  is  its  effective  coefficient  of  expansion? 

Answer: 

As  a  first  approximation,  assume  the  mass  of  each  weight  is  equal  to  its  nominal 
value.    Therefore,  the  total  volume  is  approximated  by 

.050  ,  0.030  x  0.020  _  n  ni„„  3 
WJo  +  "167F  +  TTT  ~  °-012227  cm 

The  effective  density  is  then  approximated  by 

r,  -  mm  -  .050  +  .030  +  .020     o  io  n/„m3 
P  =  M/V  *  0.012227  =  g' 

The  effective  coefficient  of  expansion  is  a  weighted  average—the  coefficients  of 
each  weight  are  added  in  proportion  to  the  volume  of  the  weight  to  the  total  volume: 

20  x  icf6  XM50  +  20  x  1Q-6  x  M30  +  69  x  10-6  x  0^20 

 16,6  16J5  2.7 

a  =   

.050  +  0.030  0.020 


16.6       16.6  2.7 


=  50  x  10"6/°C 


After  mass  values  for  the  weights  comprising  M  have  been  found,  the  metrologist 
should  verify  that  the  approximations  used  to  estimate  p  and  a  were  adequate.  Only 
very  rarely  will  this  not  be  the  case,  requiring  an  iteration:    the  mass  values 
obtained  for  the  weights  comprising  M  are  used  to  find  new  effective  values  of  p  and 
a;  these  latter  two  values  are  then  used  to  calculate  improved  mass  values. 
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4.    PROGRAM  TYPES 


We  will  discuss  two  distinct  programs  in  mass  metrology:    surveillance  and 
calibration.    As  the  name  implies,  surveillance  attempts  to  monitor  a  calibrated 
assembly  or  set  of  weights.    Surveillance  looks  for  signs  that  one  or  more  members 
of  the  set  may  have  changed  since  the  last  calibration.    Calibration,  of  course, 
attempts  to  assign  the  best  possible  value  of  mass  to  a  weight  by  relying  on  a 
chain  of  traceability  to  universally  recognized  standards.    Assignment  of  a  realis- 
tic uncertainty  to  a  calibration  value  is  equal  in  importance  to  the  calibration 
value  itself.    Several  additional  definitions  will  be  useful. 

A  standard  mass  ST  was  already  defined  in  Section  3.    We  also  have: 

SC:    Check  standard  weight  with  true  mass  (SC)T. 

This  mass  is  basically  required  to  monitor  in-house  procedures  and  accuracies. 
Usually  (SC)  is  a  calibrated  standard  S  that  has  been  added  to  a  set  during  cali- 
bration.   By  measuring  (SC)  (treated  as  an  unknown  M)  versus  S  and  comparing  the 
results  with-  the  known  value  one  can  then  determine  if  the  measuring  process  is 
either  in  or  out  of  control.    This  type  of  self  check  should  be  done  periodically 
and  the  values  for  (SC)  should  be  monitored  continuously  via  trend  charts  (see 
below) . 

Set:    A  group  of  weights  in  a  specific  order.    To  be  of  maximum  utility,  sets  usually 
cover  several  decades  of  nominal  mass;  e.g.,  1  mg  to  100  g.    Further  smaller  groups 
(called  subsets)  of  the  complete  set  can  be  selected  to  produce  any  nominal  value 
within  the  span  of  the  set  in  increments  of  the  smallest  set  members.    What  this 
means  is  that,  taking  the  example  of  a  100  g  to  1  mg  set,  any  of  the  nominal  values 
0.001  g,  0.002  g,  0.003  g,         99.998  g,  99.999  g,  100  g  can  be  obtained  by  select- 
ing appropriate  weights  from  the  set.    Usually  a  set  will  employ  the  smallest 
number  of  weights  possible  to  cover  each  decade.    A  set  with  a  sequence  5,  3,  2,  1 
could  for  instance  range  from  100  g  to  1  mg  with  the  sequence  repeated  for  each 
decade  (e.g.  500  mg,  300  mg,  200  mg,  100  mg). 

Design:    This  specifies  explicitly  how  and  in  what  order  a  set  of  weights  is  to  be 
compared  (measured).    It  states  in  detail  what  weights  of  the  set  have  to  be  used  at 
a  particular  step. 

ST:    Transfer  standard  weight  with  true  mass  (ST)T. 

As  mentioned  under  "Set"  above,  in  many  cases  a  number  of  weights,  e.g.  from 
1  kg  to  10  g  must  be  calibrated.    These  are  usually  divided  into  decade  groupings; 
that  is,  weights  from  1  kg  to  100  g  are  calibrated  at  the  same  time  using  one  design 
(A  standard  kilogram  is  also  included  in  the  design).    Next  the  weights  from  100  g 
to  10  g  are  cal ibrated—possibly  using  a  different  balance.    In  this  procedure,  it 
is  common  for  the  100  g  weight  which  was  calibrated  as  a  result  of  the  first  series 
of  measurements  to  become  the  standard  for  the  second  series.    The  100  g  weight  thus 
serves  as  a  "transfer  standard"  because  it  is  used  to  transfer  our  knowledge  of  a 
standard  kilogram  to  weights  of  smaller  denomination, 

4.1.  Surveillance 

Two  distinct  types  are  outlined  below  as  Type  I  and  Type  II.    (For  more  com- 
plete detailed  discussions  of  many  different  designs,  the  reader  is  referred  to  ref. 
[4]). 
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The  basic  idea  of  surveillance  testing  is  to  ensure  the  self -consistency  of  the 
weight  set.    For  example,  the  20  g,  30  g,  and  50  g  weights  can  be  checked  against 
the  100  g  weight  to  see  if  the  difference  is  within  expected  limits.    It  is  also 
advantageous  (though  not  essential)  to  compare  one  weight  of  the  set  (usually  the 
largest)  against  an  independent  standard,  S.    The  latter  operation  establishes 
whether  the  entire  weight  set  has  undergone  a  change--even  though  there  may  still  be 
self-consistency  within  the  set.    The  basic  motivations  for  surveillance  testing 
are: 

•  to  verify  the  values  of  newly  calibrated  weights 

•  to  establish  the  stability  of  a  new  weight  set 

•  to  determine  if  an  accident  (such  as  being  dropped  on  the  floor)  has  damaged 
the  weights  involved. 

Buoyancy  corrections  may  not  be  needed  in  surveillance  testing.  One  should 
check  the  magnitude  of  such  corrections  compared  to  the  surveillance  limits  (see 
below)  to  see  whether  it  is  worthwhile  to  make  the  buoyancy  corrections. 

4.1.1.    Type  I 

The  object  of  this  surveillance  method  is  to  perform  an  intercomparison  of  all 
weights  in  a  set  using  a  minimum  number  of  steps.    It  is  preferable  to  have  one 
standard  weight  as  a  member  of  the  set.    We  will  denote  it  as  S  for  this  discussion. 
The  standard  S  should  always  be  as  large  as  the  largest  member  of  the  set;  or  as 
large  as  is  convenient. 

For  the  comparison  measurement  one  always  starts  out  with  the  largest  weight 
and  compares  it  with  a  summation  of  weights  next  in  magnitude  such  that  the  sum  is 
equivalent  to  the  largest.    Next,  a  weight  from  the  first  summation  is  compared  with 
a  lower  summation,  and  the  process  is  continued  until  all  the  weights  have  been 
used.    If  a  standard  is  included  in  the  set,  then  the  S  has  to  be  compared  first 
with  the  largest  mass.    An  example  will  help  to  visualize  the  procedure. 

Example: 

The  set  contains  a  standard  weight  S  whose  nominal  value  is  100  g  and  weights 

of  100,  50,  30,  20,  10,  5,  3,  2,  1,  0.5,  0.3,  0.2,  0.1,  0.05,  0.03,  0.02,  0.01, 

0.005,  0.003,  0.002  and  0.001  g;  called  M-|  to  M21 ,  respectively.    Such  a  set  is 

known  as  a  set  with  mass  ratios  of  5,  3,  2,  1  from  100  g  to  1  mg.  For  the  measure- 
ment sequence,  we  start  out  with  the  standard  S  versus  M-|  =  100  g  and  work  down  to  1 

mg  such  that  all  masses  are  included  via  a  minimum  number  of  steps.    The  designation 

M.  refers  to  the  nominal  value  of  the  ith  weight  whereas  mT  refers  to  the  true  mass 

value  of  the  ith  weight. 

1st  Meas:    m{  -  ST  =  6-j ,  where  M]  =  100  g,  S  =  100  g 

2nd  Meas:    mJ  -  M1/  =  62,  where  M2'  =  (M2  +  M3  +  \) 

=  (50  +  30  +  20)  g 
=  100  g 
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T  T 

3rd  Meas:    M^  -  M^'  =  63,  where  M3'=  (Mg  +  Mg  +  M?  +  Mg) 

=  (10  +  5  +  3  +  2)  g 
=  20  g 

4th  Meas:    Mg  -  mJ'  =  64>  where  M4'  =  (Mg  +  M]0  +  M^  +  M]2) 

=  (1  +  0.5  +  0.3  +  0.2)  g 
=  2  g 

5th  Meas:    mJ2  -  mJ'  =  6g,  where  Mg1  =  (M]3  +  M]4  +  M]g  +  Mlg) 

=  (0.1  +  0.05  +  0.03  +  0.02)  g 
=  0.2  g 

6th  Meas:  m{6  -  M^1  =  6g,  where  Mg'  =  (M]?  +  M]8  +  M]g  +  M2Q) 

=  (0.01  +  0.005  +  0.003  +  0.002)  g 
=  0.02  g 

7th  Meas:    Mjg  -  mJ'  =  6y,  where  M?'  =  (M2Q  +  M2]) 

=  (0.002  +  0.001)  g 
=  0.003  g 

Note  that  the  last  measurement  is  amended  because  only  M,^  is  left  to  measure. 

With  a  very  simple  software  routine,  the  differences  (<5's)  can  then  be  compared 
against  the  known  (accepted)  values.    These  new  differences  should  then  be  plotted 
and  compared  chronologically  with  previous  tests.    Together  with  predetermined 
uncertainty  limits  one  can  then  monitor  the  particular  weight  set.    Usually  one 
devotes  one  chart  for  each  6.    For  example,  fig.  5  shows  a  survellance  chart  for  6.. 
The  horizontal  line  is  the  value  of  64  inferred  from  the  most  recent  calibration 

report  of  the  set.    The  points  represent  values  of  6^  which  were  derived  from 

surveillance  testing.    As  long  as  the  points  remain  within  the  upper  and  lower 
horizontal  lines,  (known  as  the  surveillance  limits),  we  have  no  evidence  that  any 
of  the  weights  involved  in  the  measurement  of  6^  has  changed  from  its  reported 

calibration.    If  an  obvious  trend  is  apparent,  however,  a  recal ibration  can  be  done 
before  the  surveillance  limit  is  exceeded.    Two  questions  remain  to  be  answered: 
(1)  How  are  the  surveillance  limits  determined?;  and  (2)  If  a  point  lies  outside  the 
surveillance  limits,  how  does  one  determine  which  of  the  weights  (that  is  2  g,  1  g, 
500  mg,  300  mg,  or  200  mg  in  the  case  of  fig.  4)  has  changed?    The  answers  to  these 
questions  will  be  deferred  until  after  the  description  of  Type  II  Surveillance. 

4.1.2.    Type  II 

For  Type  II  surveillance,  a  more  sophisticated  approach  to  the  mass  comparisons 
is  taken.    Again,  the  procedure  is  best  described  by  an  example. 

Suppose  we  consider  the  following  three  weight  groupings  from  the  set  used  in 
the  example  for  Type  I. 
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Weight  Grouping  Designation 
3  g  M7 

2  g  +  1  g  M3  =  M8  +  M9 

2  g  +  500  mg  +  300  mg  +  200  mg  M£'  =  Mg  +  M]0  +  M,,  +  M]2 

For  Type  II  surveillance,  the  following  three  mass  differences  are  measured. 

1.  H]-nl»  =6, 

2.  My  -  vJ3"'    =  Sz 

3.  Mj"  -  Mj" '  =  63 

(Note  that  63  cannot  be  measured  by  transposition  weighing  because  Mg  is  common  to 
M3"  and  M3,M  .) 

The  three  weighing  operations  shown  above  form  a  simple  "design".    The  concept 
of  weighing  designs  is  crucial  to  calibration  operations  and  will  be  discussed  at 
length  in  Section  4.2.2.    Here  we  may  simply  note  that  we  are  working  with  a  "3-1 ' s 
design"  (i.e.,  3  different  weight  groupings  with  1  nominal  mass). 

If  the  weighing  process  had  no  scatter  (standard  deviation  equal  to  zero),  then 
63  =  §2"         In  tnat  case  we  would  derive  no  additional  information  from  the  63 

measurement  and  its  inclusion  would,  in  fact,  be  a  waste  of  time.    Of  course,  the 
process  does  have  some  scatter  and  the  63  which  we  measure  has  only  a  slight  prob- 
ability of  being  equal  to  the  measured  value  fig  ~  61 •    Therefore,  far  from  being 

useless,  the  measurement  of  the  trio  6,,  6«>  ^3  can  be  usec'  t0  9°ocl  advantage  in  two 
ways: 

First,  a  statistical  technique  known  as  "least  squares"  fitting  [3,8]  of  the 
data  provides  better  estimates  of  the  three  mass  differences  than  can  be  gotten  with 
a  single  measurement. 

1 .  M7  -  M3"  =  1/3(26-,  +  62  -  63)  =  6-j 

2.  mI  -  M3" '  =  1/3(6-,  +  262  +  63)  =  6£ 

3.  M3"  -  M3"'  =  l/3(-6-|  +  62  +  263)  =  6^  . 

The  least-squares  estimates  6],  62,  and  63  can  be  checked  against  the  calibration 
report  as  in  Type  I  surveillance  (see  refs.  [3]  and  [8]). 

Second,  the  differences  between  the  least-squares  estimates  and  the  observa- 
tions, f<$]-<5-|  »  |V"62l  and  l63,_<53l  '  are  related  t0  scatter  ^  the  measurement 
process.6    If  one  of  these  differences  is  markedly  higher  (say  three  times  higher) 


6It  is  characteristic  of  the  least  squares  solution  to  a  3-Vs  design  that  |6-j -6 J 
=  |62-62|  =  |<$3-<$3|.    This  feature  cannot  be  generalized  to  other  designs. 
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than  the  measured  standard  deviation  of  the  balance,  one  should  redo  the  measurements 
of  6p  62,  and  63.    If  trouble  persists,  it  is  a  good  indication  that  the  precision 

of  the  balance  (eq.  (21))  has  deteriorated. 

As  another  example  we  will  use  a  set  similar  to  the  one  discussed  under  Type  I, 
i.e.  mass  ratios  5,  3,  2,  1  from  100  g  to  1  mg  with  S  =  100  g. 

The  first  3  -  1's  weighing  series  consists  of 

la.    ST  -  mJ  =  6-j ,  where      =  100  g 

lb.    ST  -  M2'  =  62,  where  M2'  =  (M2  +  M3  +  M4) 

=  (50  +  30  +  20)g 


100  g 


1c.    Mj  -  mJ'  =  63 
The  next  3  -  1's  weighing  series  is 

2a.    M3  -  M3'  =  <54,  where  M3  =  30  g 

M3'  =  (M4  +  M5) 
=  (20  +  10)g 
=  30  g 

2b.    M3  -  mJ'  =  65,  where  M4'  =  (M4  +  Mg  +  M?  +  Mg) 

=  (20  +  5  +  3  +  2)g 
=  30  g 

2c.    Mj-  -  MJ'  =  66 
This  is  followed  by  the  third  3  -  1's  series  as 
3a.    M7  -  mJ'  =  6y,  where  M?  =  3  g 

=  (2  +  1)  g 

=  3  g 

3b.        -  mJ'  =  68,  where  Mg'  =  (Mg  +  M]0  +  Mn  +  M12) 

=  (2  +  0.5  +  0.3  +  0.2)  g 

=  3  g 

3c.    mJ'  -        =  6Q 
5        6  9 
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For  the  fourth  sequence  we  have: 

4a.    m|i  -  M^'  =  61Q,  where        =  0.3  g 

V  =  <M12  +  M13> 
=  (0.2  +0.1)  g 

=  0.3  g 

4b.    Mj1  -  Mj'  =  6n,  where  Mg'  =  (M]2  +  M]4  +  M]5  +  M]6) 

=  (0.2  +  0.05  +  0.03  +  0.02)  g 
=  0.3  g 

4c.    mJ'  -  mJ'  =  6]2. 
For  the  fifth  series  we  measure: 

5a.    Mj5  -  mJ1  =  613,  where  Mlg  =  0.03  g 

M9'  =  (Mlfi  +  M1?) 

=  (0.02  +  0.01)  g 
=  0.03  g 

5b.    M{5  -  m{q'    =  6U.  where  M-j q  =  (M]6  +  M]8  +  M]g  +  M2Q) 

=  (0.02  +  0.005  +  0.003  +  0.002)  g 
=  0.03  g 

5c.    Mj'  -  M{0'  =  615. 
Finally,  we  arrive  at  the  last  series  with 
6a.    m|9  -  mL  =  616  where  M]9  =  0.003  g 

Mll*  =  (M20+M21> 

=  (0.002  +  0.001)  g  =  0.003  g 

6b.    Mjg  -  Mj2'  =  617  where  M^'    =  (M2Q  +  M22) 


=  (0.002  +  0.001*)  g  =  0.003  g 

u18' 


6c.    Nil'  -  mIa  =  6 


Note  that  in  the  last  series  a  known  (standard)  0.001-g  weight  was  added  in  order 
to  complete  the  design. 
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Once  again  the  least  squares  fitting  technique  provides  better  estimates  such 
that  eqs.  (la),  (lb),  and  (lc)  yield 

la')    1/3(261  +  <$2"63)  =  61 ' 
2b')    1/3(61  +  262  +  63)  =  6£ 
3c')    l/3(-61  +  62  +  2  63)  =  63. 

If  for  any  reason  a  suitable  standard,  S,  is  unavailable,  then  any  uncalibrated 
weight  or  summation  of  weights  of  the  correct  nominal  value  can  be  used  in  eqs.  (la) 
and  (lb).  However,  in  this  case  only  6i  can  be  checked  against  the  reported  (known) 
calibration  and  only  internal  consistency  of  the  entire  set  can  be  determined. 

We  also  have 


2a') 

l/3(264  +  65  - 

«6)  ■  V 

2b') 

1/3(6-  +  26,  + 

O  =  V 

0  0 

2c') 

and 

3a') 

l/3(-64  +  65  + 
l/3(26?  +  6g  - 

266>  =  66 
69)  =  6) 

3b') 

l/3(6y  +  26g  + 

6g)    =  6g 

3c') 

l/3(-67  +  6p  + 
/  0 

26q)  =  6' 

as  well  as 

4a') 

l/3(26ln  +  6,, 

1  U  II 

-  <512)  = 

4b') 

l/3(6ln  +  26n 

+  619)  = 
1  c 

1  1 

4c') 

l/3(-6,0  + 

+  2«12)  = 

612 

Finally 

5a' ) 

l/3(2«,3  +  «M 

613 

5b') 

1/3(613  +  2614 

+  «15>  = 

514 

5c') 

1/3(-613  +  <514 

♦  2«,5)  = 

fi15 

and 

6a') 

1/3(2«16  +  617 

"  618'  = 

516 

6b') 

V3(S16  +  261? 

+  618>  * 

617 

6c') 

l/3(-616  +  61? 

+  2«18)  = 

618 

The  differences  obtained  for  the  first  two  results  of  each  3-1 's  series  (i.e.  a 
and  b)  should  be  compared  to  accepted  values  derived  from  the  most  recent  calibration 

report. 
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The  above  examples  of  surveillance  measurements  for  Type  I  and  Type  II  present 
a  simplified  picture  of  the  equations.    All  6's  must  contain  any  necessary  buoyancy 
corrections  for  true  mass  differences.    Different  6's  are  obtained  if  apparent 
masses  are  used.    From  the  simplified  eq.  (19a),  one  finds 

MA  -  SA  -  kA6  +  (pA  -  pQ)  (VQM  -  Vos)  =  6 

where  we  assumed  that  the  mass  of  the  standard  is  given  in  apparent  mass  units,  S  . 

Deviations  of  the  measurements  from  values  inferred  from  the  most  recent  cali- 
bration report  are  evident  from  the  previously  mentioned  surveillance  charts. 

4.1.3.    Surveillance  Limits 

Upon  carrying  out  a  surveillance  test  of  either  Type  I  or  II  it  will  become 
evident  that,  after  any  necessary  buoyancy  corrections  have  been  made,  the  measured 

values  6.  (Type  I)  or  the  predicted  values  6!  (Type  II)  do  not  exactly  agree  with 

those  values,  fi.9  calculated  from  the  calibration  certificate  accompanying  the  set. 

To  judge  how  serious  the  disagreement  is,  one  must  also  calculate  surveillance 
limits.    The  surveillance  limits  associated  with  each  6..  or  61  define  the  limits  of 

credibility  that  the  difference   6^  -  6.    or   6^  -  6.j    could  be  due  to  a  combination 

of  calibration  uncertainties  of  the  weight  set  and  the  random  error  of  the  measure- 
ment which  is  associated  with  the  balance  used  in  the  surveillance  measurements. 

Let  us  designate  to  random  and  systematic  error  limits  by  SL; 

(SL)  =  U  +  E  (42) 

where 

U  =  systematic  error  as  determined  from  the  calibration  report 
E  =  limit  to  random  error  =  3a 
and  o  =  standard  deviation  of  the  measurement  (refer  to  Appendix  D  for 

information  on  estimating  this  number). 

Then  any  value  6.  or  6!  should  fall  between  6?±SL,  as  shown  in  fig.  5.    The  value  of 
3a  gives  a  99.7  percent  confidence  level  for  the  random  error  (i.e.  if  the  measure- 
ments were  repeated  a  great  many  times  there  is  a  99.7  percent  chance  the  average 
would  be  within  ±  E  of  the  result  of  a  single  surveillance  test). 

For  surveillance,  U  is  the  root-sum-square  of  the  individually  reported  cali- 
bration uncertainties  such  that 

U  =  (EU?)1/2  (43) 

Combining  the  uncertainties  in  this  way  is  strictly  valid  only  if  they  are  uncorre- 
cted.   It  is  the  nature  of  calibration  designs,  however,  that  the  uncertainties  of 
weights  within  the  set  generally  are  correlated.    Thus  one  should  view  eq.  (43)  as  an 
approximation  which  is  adequate  for  surveillance  limits.    The  reader  is  referred  to 
the  discussion  of  calibration  uncertainties  below  for  a  fuller  explanation. 

The  standard  deviation  used  is  known  from  many  previous  measurements  (eq.(21)). 
It  is  a  measure  of  the  random  errors  in  the  balance  being  used. 

At  this  stage  it  is  worthwhile  to  quote  a  couple  of  examples  and  step  through 
the  process.    From  Type  I  surveillance,  we  have 
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S  =  10  g 
M,  =  100  g 

=  100  g  =  50  g  +  30  g  +  20  g  =  M2  +  M3  +  M4 
Assume  the  following  previously  reported  calibration: 

Mass  Value  Uncertainty 
ST  =  100.0010196  g  Us  =  0.015  mg 

NlJ  =  100.0009407  g  UM  =  0.020  mg 

M2  =    50.0004628  g  =  0.011  mg 

M3  =    30.0002926  g  =  0.012  mg 


NlJ  =    20.0001578  g  UM    =  0.010  mg 


Also  assume  that  a  =  0.026  mg. 
Since  the  first  measurement  is 


M 


j  -  S1  =  6p  6^  =  -0.0000789  g  =  -0.079  mg 


and  we  have  t 

U  =  \U 


2  +  U  2 
S  UM1 


=  \(0.015)2  +  (0.020)2  =  0.025  mg 
and  (SL)  =  U  +  3a 

=  0.025  +  3  x  0.026 
=  0.103  mg. 

Hence  (mJ  -  ST)  values  should  fall  in  the  range  of  (-0.079  ±  0.103) 
For  the  second  measurement  we  have 


so  that 


m{  -  M^'  =  62;  62  =  0.028  mg 


Jm.  =  Jul  +  U?,  +  U? 

M£      yj  K,  M3 


=  0.019  mg 


=  ^ (0.019)2  +  (0.020)2 


=  0.028  mg 
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and  (SL)  =  0.028  +  3x0.026 

=  0.106  mg 

and  mI  -  mJ'  should  fall  in  the  range 

0.028  ±  0.106  mg. 

Next,  we  consider  the  example  from  Type  II.  The  first  series  of  the  design  consists 
of  S,  M-j ,  and  M£  with  M£  =  Nl>  +  M3  +  M^.    Let  us  assume  the  same  calibration  report 

as  above,  with  the  scale  standard  deviation  at  a  =  0.026  mg.    In  this  case  we  have 
three  measurements 

ST  -  m|    =  6] 

sT  -  mJ'  =  <$2 
m{  -  mJ-  =  63 

Because  we  have  used  a  3  -  Ts  design,  the  standard  deviation  of  our  result  is  not 
simply  the  standard  deviation  of  a  single  measurement,  a,  but  is  instead  }/2/Z  o.  A 
proof  of  this  conclusion  is  beyond  the  scope  of  this  text  but  is  explained  fully  in 
ref.  [3]. 

There  are  several  surveillance  limits  that  can  be  calculated.    For  the  three 
differences,  we  follow  the  calculations  in  Type  I  (but  we  have  to  remember  that  we 
are  in  Type  II  so  that  the  standard  deviation  is  /2/3  a). 

Limits  for  ST  -  mI 


U  = 


We  then  have 


V  us  +  UM    =^(0.01^  +(o.02<^    =  0.025  mg 


so  that 


Limits  for  ST  -  M^' 


(SL)  =  U  +  3  x  a)  =  0.025  +  3/273  x  0.026 

=  0.089  mg 

(ST  -  mJ)  =  (0.079  ±  0.089)  mg. 


where  UM'  =   V  U*j   +^    +  U|/    =  V^0-011)2  +  (°-012)2  +  (0.010): 


so  that 
Hence 


1  "2  "3 
0.019  mg 


U  =    y (0.015)2  +  (0,019)2    =  0.024  mg 


(SL)  =  U  +  3/2/3  o  =  0.024  +  3/273  0.026  =  0.088  mg 
and  (ST  -  mJ')  =  (0.106  ±  0.088)  mg. 
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4.1.4.    Identifying  Weights  Which  Have  Changed 

If  a  measurement  falls  outside  the  surveillance  limits,  it  becomes  necessary  to 
determine  which  individual  weights  are  responsible.    Simple  deductive  reasoning  is 
all  that  is  required,  although  a  few  extra  weighing  combinations  may  also  be  needed 
if  Type  I  surveillance  was  used. 

As  an  example,  suppose  that  in  the  Type  I  measurements  shown  earlier  (4.1.1)  the 
value  for  64  was  outside  but  the  value  for  63  was  inside  the  surveillance  limits. 

All  we  know  at  this  point  is  that  there  is  some  change  in  the  subset      =  1  g  +  0.5  g 

+  0.3  g  +  0.2  g.    Three  additional  weighings  are  now  required  to  pinpoint  the  cause 
of  the  discrepancy.    Recall  that  MQ  =  1  g,  Mln  =  0.5  g,  M-,,  =  0.3  g,  M19  =  0.2  g, 
M13  =  0.1  g.  y  lu  11  u 

We  then  measure 

4'  )    mJ  -  Mj"  =  6J,  where  MJJ  =  M]0  +        +  M]2 

4"  )    Mj0  -  M5"  =  65,  where  M£  =  Mn  +  M]2 


4'")          -  Mg"  =  6£,  where  M£  =  M]2  +  M 


13 

After  the  surveillance  limits  for  6j|,  6£,  6£  are  calculated  one  investigates  the 
results: 

1.  )    If  6^  lies  outside  the  surveillance  limits  but  6£  and  6£  do  not,  then  it  is 

probable  that  Mg  has  changed  since  its  last  calibration. 

2.  )    If  <5^  and  6£  lie  outside  the  surveillance  limits  by  opposite  amounts  and  6£  is 

inside,  then  it  is  probable  that  the  M-|q  (0.5  g)  weight  has  changed. 

3.  )    If  6^  and  6£  lie  outside  the  limits  by  about  the  same  amount  and  6£  is  outside 

also  by  the  same  amount  but  in  the  opposite  direction,  then  it  is  probable  that 
the        (0.3  g)  weight  has  changed. 

4.  )    If  64,  6£,  and  6£  all  lie  outside  the  limits  by  about  the  same  amount  in  the 

same  direction,  then  it  is  probable  the  M-j2  (0.2  g)  has  changed. 

If  none  of  the  above  conditions  are  met,  then  it  is  probable  that  more  than  one 
weight  has  changed.    The  reader  should  then  consult  ref.  [4]  for  a  more  thorough 
analysis  of  surveillance  methods. 

4.2.  Calibration 

The  process  of  "calibration"  assigns  mass  values  to  weights  by  comparing  the 
unknown  weights  to  recognized  standards.    An  uncertainty  1 imi t — 3  times  the  standard 
deviation  of  the  measurement  process,  plus  estimated  uncertainties  systematic  to  the 
measurement  process—accompanies  each  calibrated  value.    To  an  even  greater  degree 
than  in  Type  II  surveillance,  redundant  information  is  gathered  in  order  to  determine 
whether  the  measurement  scatter  is  acceptable.    A  powerful  self-consistency  check  of 
the  calibration  process  is  also  included. 
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Most  of  the  following  information  can  be  found  in  much  greater  detail  in  refs. 
[1]  and  [3].  It  is  repeated  here  mainly  for  explaining  the  overall  approach  of  the 
program. 

Weight  sets  come  in  various  denominations.    The  following  groupings  are  fre- 
quently calibrated: 

1.    Nominally  Equal  Groups 

These  sets  consist  of  weights,  all  of  which  have  the  same  nominal  value.  The 
number  of  weights  can  go  from  3  to  as  high  as  13  or  more.    (Sets  of  up  to  50 
members  have  been  utilized.)    For  large  sets,  subsets  are  usually  analyzed. 
Otherwise  the  number  of  weighing  operations  would  become  impractical ly  large. 


2.    Groups  for  2,  2,  ...  1,  1 


sets. 


Several  combinations  are  in  use  such  as  2,  1,  1,  1,  or  2,  2,  1,  1,  1,  or  2,  1, 
1,  1,  or  2,  2,  1,  1,  1,  etc.    Here  2  and  1  imply  the  same  decade  of  mass  such 
as  2  kg  and  1  kg. 

3.  Binary  and  Miscellaneous  Groups 

These  sets  are  usually  comprised  of  pound  units.    Many  combinations  are  available, 
the  use  of  which  is  slowly  diminishing.    A  few  examples  provided  in  ref.  [3] 
are  4,  3,  2,  1,  1,  or  10,  5,  2,  2,  1,  1,  or  6,  5,  4,  3,  2,  1,  etc.    Note  that  a 
combination  of  the  lower  masses  always  sums  to  the  value  of  a  higher  mass  unit. 

4.  The  5,  3,  2,  1  and  5,  2,  2,  1  Groups 

These  are  the  sets  most  widely  used.    In  general  a  set  of  weights  in  this  group 
can  span  many  decades  with  each  decade  comprising,  for  example,  the  5,  3,  2,  1 
sequence.    In  many  cases  standards  or  "check  standards"  (see  below)  are  added 
to  the  sets.    Examples  of  such  groups  could  be  5,  5,  3,  2,  1,  1,  or  5,  3,  2,  1, 
1 ,  etc. 

4.2.1.    Trend  Elimination  in  Direct  Reading  Balances 

If  a  single-pan,  direct  reading  balance  has  truly  constant  sensitivity  and  is 
subject  to  only  a  slight  linear  drift,  the  following  simplification  is  possible: 

Suppose  there  are  four  nominally  equal  weights,  designated  A,  B,  C,  and  D  to  be 
compared  in  a  particular  calibration  scheme,  the  following  eight  "direct"  weighings 
are  done: 


(1 
(2 
(3 
(4 
(5 
(6 
(7 
(8 
(9 
(10 

(11 
02 


Place  A  on  balance  and  read  0-| 
Remove  A 

Place  B  on  balance  and  read 
Remove  B 

Place  C  on  balance  and  read 
Remove  C 

Place  D  on  balance  and  read  6, 
Remove  D 

Replace  D  on  balance  and  read  65 
Remove  D 

Place  C  on  balance  and  read  6^ 
Remove  C 


(13)    Place  B  on  balance  and  read  6^ 
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(14)  Remove  B 

(15)  Place  A  on  balance  and  read  6g. 

We  have  taken  the  weights  in  the  following  order:    A,B,C,D,D,C,B,A.    We  then 
make  the  following  calculations: 

eA  =  1/2(0-!  +  e8) 
eB  =  i/2(e2  +  e7) 
ec  =  i/2(e3  +  e6) 
eD  =  i/2(e4  +  e5) 

Then  the  estimated  difference  in  mass  between  any  two  weights  (say  C  and  D)  is 

CT  -  DT  =  pA(vc  -  vD)  +  k(ec  -  eD) 

where  we  assume  k  is  well-known  and  constant.  The  symmetry  of  the  weighing  sequence 
removes  problems  caused  by  any  linear  drifts  in  the  balance  or  weighing  conditions. 

The  above  example  is  an  instance  of  "trend  elimination."    (Double  substitution 
is  another  instance.)    When  weighings  are  subject  to  large  air  buoyancy  corrections, 
drifts  in  temperature  or  barometric  pressure  may  lead  to  errors  if  those  quantities 
are  only  read  once  during  all  the  measurements  required  by  the  design.    Some  cali- 
bration schemes,  however,  have  the  property  of  "trend  elimination"  even  for  this 
problem  [3]. 

4.2.2.  Designs 

Sets  or  subsets  of  weights  are  calibrated  together  by  means  of  a  weighing  de- 
sign.    A  design  simply  prescribes  what  weighings  are  to  be  made.    Each  weighing  is 
used  to  estimate  a  mass  difference  between  two  nominally  equal  weights  or  groups  of 
weights  in  the  set. 

A  typical  design  for  a  group,  say  5,  3,  2,  1,  1,  1,  could  be  as  follows: 

Mass 

Observation  M5     M3     M2     M1     M£1  SC 

61  +      -       -  + 

62  +  + 

63  +       -       -  - 

64  +       -  - 

65  +  .... 

66  +       -       +       -  - 

67  +       -       -  + 

68  +       -       -       -  + 
6g  *       -  - 

610  +  " 

611  "  " 
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Here  +/-  signs  indicate  the  weight  starts  out  on  the  left/right  pan,  for  transposi- 
tion weighing,  or  that  the  weight  is  the  first/second  in  substitution  weighing.  We 
notice  that  in  this  group  there  are  six  masses,  which  must  be  determined  (K  =  6),  11 
measurements  (N  =  11),  and  suppose  1  restraint  (e.g.  the  mass  of  the  5+3+2  summation 
is  known  from  previous  measurements,  L  =  1).    We  can  then  calculate  the  number  of 
degrees  of  freedom  (D.F.)  for  this  design  as 

D.F.  =  N-K+L  =  11-6  +  1  =  fr. 

From  the  11  observations  and  the  given  value  of  the  restraint  one  can  then  use  the 
least  squares  method  to  solve  all  equations  to  obtain  best  values  of  all  differences 
and  therefore  all  masses,  standard  deviations,  and  variances.    Such  fitting  techniques 
are  well  described  in  literature  available  from  the  NBS  [3].    (See  Appendix  D.) 

In  general  the  more  degrees  of  freedom  provided  by  a  particular  design,  the 
lower  the  estimates  of  standard  deviation.    Recall  the  previous  section  on  surveil- 
lance testing.    Each  measurement  in  Type  I  surveillance  had  no  degrees  of  freedom. 
That  is,  we  needed  to  determine  a  quantity  6.  and  we  had  only  one  measurement  of  that 
quantity  (K  =  1 ,  N  =  1 ,  L  =  0;  D. F.  =  0) .    In  Type  II  surveillance,  however,  we  used 
a  simple  3-1  design:    Three  quantities  were  determined,  for  example  6-j,  6^,  $3 

(N  =  3),  three  measurements  were  made  (6-j,        63 »  so  that  K  =  3),  and  there  was  one 

restraint  (6i  =  6£  -  6-j).    Thus  D.F.  =3-3  +  1=1.    The  appearance  of  a  degree  of 

freedom  allowed  a  crude  statistical  check  (e.g.,  whether   6^  -  6-j    <   balance  standard 

deviation)  and  also  led  to  a  slightly  reduced  standard  deviation  (/2/3  a  instead  of 
a). 

The  design  illustrated  in  this  section  is  much  more  sophisticated  than  the  3-1 
but  is  motivated  by  the  same  twin  desires  for  good  statistical  analysis  and  small 
standard  deviation  of  the  calibrated  masses.    The  design  in  the  above  example  is 
shown  in  figure  6,  where  it  is  designated  C.2.    This  figure  is  reproduced  from  ref. 
[3].    Note  that  two  different  restraints  are  considered:    restraint  A  is  that  the 

7 

mass  of  the  summation  of  the  5,  3,  and  2  weights  is  known;    restraint  B  is  that  the 

o 

mass  of  a  single  "1"  weight  is  the  known  standard.     We  will  discuss  the  results 
based  on  restraint  A  in  detail.    That  is,  we  will  assume  that  the  true  mass  of 

M5+M3+M2  is  known*    Let  us  designate  this  mass  as  rT- 

In  the  532111  design  shown  above,  we  have  labeled  one  of  the  weights  as  Zl . 
This  indicates  that  we  could  calibrate  the  5+3+2  summation  of  the  next  lower  decade 
of  the  same  weight  set  which  could  then  serve  as  a  transfer  standard  for  a  subsequent 
design.    By  "SC"  we  designate  a  weight  which  is  external  to  the  set  being  calibrated 
and  which  will  serve  as  check  standard. 


This  would  be  useful,  for  instance,  if  the  5+3+2  summation  had  been  calibrated  in  a 
previous  design.    In  this  case,  the  5+3+2  summation  would  be  a  transfer  standard. 
This  restraint  is  used  in  working  down  in  mass  from  1  kg. 

8This  would  be  useful  in  working  up  in  mass  from  1  kg. 
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The  least-squares  solutions  for  the  five  unknown  weight  are  [3]: 
=  (1/920)  {100(61+62+63+64)+6065  -  20(66+67+68+6g+610+61  ]  )+460RT} 
M3  =  (l/920){-68(61+62+63+64)-465+124(66+67+68)-60(6g+610+611)+276RT} 

M2  =  (l/920){-32(61+62+63+64)-5665-104(66+67+6g)+80(69+610+611)+184RT} 

m{  =  (l/920){11961+462-ni63+464-10865+12866-102(67+68)-125(69+610)-10611+92RT} 

=  (l/920){-11161+11962+4(63+64)-10865-10266+12867-10268-1256g-10610-125611+92RT} 
In  addition,  the  least-squares  solution  provides  a  value  for  the  check  standard: 

(SC)T  =  (l/920){461-11162+11963+464-10865-102(66+67)+12868-1069-125(610+611)+92RT} 

The  metrologist  must  determine  the  values  6.,,  62,  ...  6^  and  the  value  of  R*^  by 

experiment.    The  least-squares  solution  shows  how  to  combine  the  observations  with 
the  value  of  the  restraint  in  order  to  arrive  at  mass  values  for  the  weights  used  in 
the  design. 

In  addition  to  mass  values,  the  least-squares  solution  provides  fitted  values 
for  the  observations >  which  when  compared  with  the  measured  values  give: 


V 

&\  - 

61-M5+M3+M2-H1+ME1 

62-Hl+Mj+Mj-My1T+(SC) T 

VM5+M3+M2+M1 -(SC) 

V 

64"M5+M3+M2 

V 

■65  = 

fig-Mg+Mg+M^+M^  +  CSC) 

V 

■&6  = 

66-M3l+M2-M1l+Mzl1  +  (SC)  1 

V 

•57  = 

67-M3+M2+M1-Mzl+(SC) 

V 

•«8  = 

68-M3+M2+M1+MEl-(SC) 

V 

aJ  aT 
69-M2+M1+Mzl 

610' 

'6io= 

610-mJ+m{+(sc)t 

61T 

-«ir 

~T  ~T     ^  T 
611-M2+MZ1  +  (SC) 
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These  deviations  are  useful  for  two  reasons.  First,  the  estimated  standard  deviation 
of  the  least-squares  fit  is  equal  to 

.  1/2 

'i    (6,-6;  )2 

vi=l  11 

Second,  a  glance  at  each  of  the  eleven  values  of  6.-61  can  often  pinpoint  the  source 

of  a  blunder  in  entering  the  raw  data  into  a  computer— that  is,  if  a  blunder  has 
been  made,  the  value  of  6.-61  affected  by  the  mistake  will  often  appear  much  larger 
than  all  the  other  values. 

Finally,  the  least-squares  solution  also  provides  estimates  of  the  standard 
deviation  for  each  computed  mass  and  each  combination  of  mass.    Below  the  design 
in  fig.  6  is  a  table  called  "Factors  for  Computing  Standard  Deviations."  Suppose 
the  standard  deviation  of  the  measurement  process  is  a.    Then  if  design  C.2  is  used 
to  calibrate  the  mass  of  the  "5"  weight,  for  example,  the  table  tells  us  that  the 
standard  deviation  assigned  to  this  mass,  o^,  will  be: 

0.2331  x  a  if  restraint  A  were  used 
1.7846  x  a  if  restraint  B  were  used. 

Also  from  the  table  we  see  that  the  standard  deviation  of  the  sum  of  the  5  and  3 
weights,  Og,  is  given  by: 

0.2638  x  a  if  restraint  A  were  used 
2.8284  x  a  if  restraint  B  were  used. 

It  is  typical  of  least  squares  results  that,  for  example, 

2        2  1/2 

Tn  +l.  am+n  ^  ^am   +  an  ^ 

In  this  case, 

2        2  1/2 
°8  <  ^a5    +  a3  '       ^or  restrdint  A 

2        2  1/2 
a8  >  (°k    +  a<?  )       f°r  restraint  B. 


Thus  the  metrologist  must  choose  both  design  and  restraint  carefully  to  minimize  the 
resulting  standard  deviations. 

4.2.3.    Statistical  Checks 

Computer  programs,  such  as  those  developed  by  NBS,  are  routinely  utilized  in 
mass  calibration  laboratories.    The  user  supplies  all  the  measured  data  for  the  set, 
all  environmental  conditions,  and  other  necessary  data.    The  program  then  provides: 

1.  a  detailed  listing  of  data  provided 

2.  the  least  squares  fit,  i.e.  the  desired  mass  values 

3.  the  "F  ratio"  and  the  "t  value". 

Two  crucial  assumptions  underlie  the  calibration  of  unknown  weights  by  least- 
squares  fitting  of  design  data: 
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(Reproduced  from  ref.  [3]). 
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1.  The  scatter  in  the  data  just  taken  is  typical  of  the  scatter  found  in  pre- 
vious measurements. 

2.  The  mass  of  the  standard  weight  used  in  the  calibration  design  has  not 
changed  from  its  accepted  value. 

The  F  ratio  and  t  value  provide  important  tests  of  these  two  assumptions. 

The  F  ratio  essentially  monitors  the  precision  of  the  measuring  process.  A 
detailed  discussion  is  presented  in  Appendix  D.    We  note  that 


where  a  is  the  variance  of  the  particular  balance  being  utilized  (based  on  a  large 
collection  of  previous  measurements).  Generally,  we  expect  F  to  be  close  to  1,  but 
we  must  not  be  surprised  if  a  particular  value  of  F  is  somewhat  larger  than  1. 

In  particular,  a  simple  check  for  F  can  be  (and  has  been)  established  for  most 
mass  metrology  laboratories.    By  comparing  F  with  a  fixed  ratio  F  ,  which  could  be 

defined  at  the  99  percent  confidence  level  one  can  then  easily  monitor  whether  the 
measured  ratio  F  is  greater  than  F. 


If  this  check  holds  true,  then  the  measurement  process  is  considered  "out  of  control" 
and  further  studies  have  to  be  conducted.    The  quantity  Ft  depends  only  on  the 

2 

degrees  of  freedom  (k-j)  in  s    and  on  the  probability  level  at  which  we  wish  to  con- 
duct the  test.    The  calculations  are  outlined  in  Appendix  D. 

The  second  test  is  called  the  "t  test"  which  monitors  the  accuracy  of  the  mea- 
surements.   For  this  test,  an  extra  mass,  the  check  standard,  is  included  in  the 
particular  design.    A  simple  check  is  then  performed  to  see  whether  the  mass  value 
assigned  to  the  check  standard  by  the  calibration  agrees  with  the  accepted  value. 
The  t  test  thus  monitors  the  systematic  errors  of  the  measuring  process. 

In  particular  one  calculates 


i  .e. 


tc  =  rsc'   "  accepted  (sc)  |/ac 


where 


(SC)T  =  observed  mass  of  the  check  standard  as  found  by  least-squares  fitting; 
accepted  (SC)T  =  accepted  mass  of  the  check  standard; 


and 


lcZ  +  oT2  +  (Nsc/NR)2aR2 
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Here 

Oj  =  accepted  variance  of  the  measurement  process  between  runs  (see  below) 

I  =  multiplier  which  is  determined  by  the  least  squares  process, 

c  =  "process"  standard  deviation  of  the  balance 

N  =  the  nominal  value  of  the  check  standard 

NR  =  the  nominal  value  of  the  standard  or  transfer  standard  used  asthe 
restraint  in  the  least  squares  solution 

aR  =  the  standard  deviation  of  the  value  of  the  restraint,  oR  is  taken  as 
zero  unless  the  restraint  is  a  transfer  standard,  whose  value  was  de- 
termined as  part  of  the  calibration  of  the  complete  weight  set  [1]. 

For  simplicity,  we  will  discuss  a  case  for  which  aR  =  0.  In  general  one  checks  to 
see  that 


The  probability  is  less  than  one  in  100  that  this  inequality  will  be  violated  by 
chance.    Thus  t   =  3  is  taken  as  the  control  limit.    Any  measurement  for  which  t   >  3 
is  considered  "Out  of  control"  and  is  repeated. 

CAUTION:    The  t  test  is  one  of  the  best  statistical  measures  of  systematic 
errors  available.    Nevertheless,  by  their  very  nature,  systematic  errors  are  diffi- 
cult to  detect.    If,  for  instance,  buoyancy  corrections  were  not  important  in  assign- 
ment of  mass  to  the  check  standard  but  were  important  in  assigning  mass  values  to 
the  weights  being  calibrated,  the  t  test  will  not  detect  systematic  errors  in  the 
buoyancy  correction.    Also,  the  t  test  cannot  detect  an  identical  change  in  the 
standard  and  check  standard.    We  have  taken  a  very  simple  example  of  a  t  test.  The 
most  general  case  is  shown  in  ref.  [1]. 

2 

The  variance  of  the  measurement  process  between  runs,  Oy  ,  is  an  important 

concept  [1,8].    Recall  that  the  process  standard  deviation,  a,  is  a  measure  of 
scatter  over  the  period  of  time  it  takes  to  make  a  mass  measurement:    i.e.,  a  few 
minutes  to  a  few  hours.    There  may,  however,  be  sources  of  random  uncertainty  which 
fluctuate  more  slowly  (i.e.,  over  days  or  months);  but  still  rapidly  compared  with 
intervals  between  recal ibration  of  a  given  weight.    How  can  we  estimate  the  process 
standard  deviation  for  this  longer  period  of  time?    The  easiest  way  is  to  monitor 
the  values  of  mass  assigned  to  the  check  standard  each  time  it  is  used  in  a  run. 
After  many  runs,  over  a  span  of  many  months,  we  can  estimate  the  variance,  o  S  of 
the  observed  values  of  the  check  standard  about  their  mean  value. 

If  the  process  standard  deviation  has  no  "between-run"  component,  then 

o\  =  £o2    (o2  *  0) 

2 

where  a    is  the  process  variance  which  was  found  by  pooling  variances  estimated  from 

the  least  squares  fits  to  many  individual  runs.    That  is,  o  is  the  within-run  pro- 

2 

cess  variance  based  on  combining  calculations  of  s    from  many  runs. 
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It  will  generally  be  true,  however,  that 

2     „  2 
oc  >  to 

indicating  that  there  is  an  additional  source  of  scatter  from  run  to  run.    We  define 
this  between-run  variance  for  the  measurement  process  as 

4  =  °l  '  £°2  • 

2 

If  Oj  is  not  negligible,  the  source  of  between-run  scatter  should,  of  course, 

be  sought  and  eliminated  if  possible?  Failing  this,  the  process  standard  deviation 
must  be  expanded  so  that  the  over-all  uncertainty  assigned  in  the  calibration  report 
includes  the  between-run  component. 

The  cautionary  statement  given  above  for  the  t  test  also  applies  to  the  use  of 
the  check  standard  in  estimating  between-run  scatter.    That  is,  if  a  weight  being 
calibrated  differs  markedly  in  its  construction  from  the  weight(s)  used  in  the  check 
standard,  a  cry  component  may  go  undetected  if  one  relies  only  on  check-standard 
behavior. 

Suppose  a  series  of  design  data  is  "in  control."     True  mass  values  have  been 
assigned  to  the  unknown  weights  and  uncertainty  limits  have  been  given.    What  do 
these  uncertainties  mean?    In  a  practical  sense,  the  meaning  is  that  the  uncertainty 
bands  assigned  to  each  mass  calibrated  by  the  methods  of  ref.  [1]  should  almost 
always  overlap  the  uncertainty  bands  which  would  have  been  obtained  had  the  masses 
been  calibrated  at  NBS. 

Note  that  this  specification,  while  satisfying  the  expectation  of  the  overwhelming 
majority  of  users,  is  somewhat  less  stringent  than  asserting  that:    the  true  mass  as 
defined  in  eq.  [1]  and  reported  in  S.I.  units  is  almost  surely  within  the  uncertainty 
band  provided  by  the  calibration. 


9aT  has  been  found  to  be  negligible  for  the  NBS  calibration  process. 
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5.  CONCLUSION 


This  document  has  been  less  of  a  complete  treatise  on  mass  metrology  than  a 
guide  for  the  mass  metrologist.    We  have  indicated  appropriate  references  where  the 
various  topics  touched  on  are  treated  in  greater  detail. 

What  we  have  tried  to  offer  is  a  coherent  development  of  the  fundamental  con- 
cepts of  mass  metrology  as  a  means  of  explaining  the  relevance  to  the  metrologist  of 
the  references  found  at  the  end  of  this  publication. 

In  brief,  starting  with  the  basic  law  of  classical  mechanics,  we  derived: 

Relations  for  the  true  mass  differences  between  nominally  equal  weights 
(eqs.  (10')  and  (19b)) 

Relations  for  apparent  mass  differences  between  nominally  equal  weignts 
(eqs  (18)  and  (19a)) 

We  stressed  the  importance  of  measuring  or  calculating  air  density  in  order  to  make 
proper  buoyancy  corrections. 

We  then  provided: 

An  explanation  of  how  the  basic  mass  relations  apply  to  measurements  on 
several  types  of  commonly  used  balances. 

A  demonstration  of  how  these  relations  are  used  in  two  types  of  surveil- 
lance testing. 

0      A  brief  sketch  of  how  these  relations  are  used  in  mass  calibrations. 

We  concluded  with  a  discussion  of  the  statistical  checks  crucial  to  a  calibration. 

In  a  document  such  as  this,  there  is  a  danger  of  losing  the  basic  outlines  of  a 
mass  measurement  program  in  the  many  details  essential  to  carrying  out  a  successful 
program.    Every  program  must  have  a  goal --for  example,  the  calibration  of  weights  in 
a  certain  range  of  masses  to  a  desired  uncertainty. 

One  must  then  acquire  balances  equal  to  the  task  and  house  them  in  a  room  which 
will  not  degrade  their  performance.    Control  charts  should  be  used  to  establish  the 
longterm  reliability  of  the  balances.    Calibrated  sets  of  weights  to  be  used  as 
working  standards  should  be  acquired.    Surveillance  testing  can  establish  the  stabil- 
ity of  these  weights  to  within  the  surveillance  limits. 

Anticipated  levels  of  buoyancy  correction  should  be  estimated.    Auxiliary  equip- 
ment such  as  barometers,  hygrometers, and  thermometers  should  be  acquired,  if  necess- 
ary, to  achieve  adequate  capability  in  determining  the  density  of  air.    These  instru- 
ments must  be  calibrated  periodically. 

When  calibrating  weight  sets,  weighing  designs  should  be  used.    These  are 
selected  on  the  basis  of  providing  the  necessary  accuracy  with  the  least  number  of 
weighing  operations.    Great  care  must  be  exercised  in  assigning  a  total  uncertainty 
to  calibration  results. 

The  techniques  outlined  in  this  document  form  the  core  of  such  a  mass  program. 
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APPENDIX  A 


APPARENT  MASS  OF  BUILT-IN  BALANCE  WEIGHTS 

Balance  manufacturers  almost  always  adjust  the  built-in  dial  weights  of  single- 
pan  balances  so  that  their  apparent  masses  equal  the  value  of  the  dial  position. 
No  matter  of  what  density  material  the  manufacturer  has  actually  made  his  weights, 
the  user  generally  assumes  that  the  true  masses  equal  the  dial  values  and  the 
densities  of  the  built-in  weights  equal  pR,  the  basis  density. 

Let  us  see  why  this  scheme  works  and  what  its  limitations  are.    These  consider- 
ations can  be  demonstrated  by  referring  to  a  single  dial -weight  (of  density  pd) 
that  exactly  balances  a  weight  X,  which  is  on  the  pan. 

We  begin,  since  we  are  discussing  a  balance  of  the  type  shown  in  fig.  3  with 
eq.  (22): 

k6l  '=  (MX'PAVX)  "  (£dT-pazv.)  +  ke0, 

where  we  have  renamed  M  by      and      by  V^.    We  will  assume  that  6q  was  adjusted  to 

zero  before  X  was  placed  on  the  balance.    We  have  chosen  a  special  case:    after  X 
was  placed  on  the  balance,  the  screen  reading  6-j  returned  to  zero  with  the  removal 

of  a  single  dial  weight.    Let  us  refer  to  this  dial  weight  as  "D"  and  say  that  it 

has  a  true  mass      and  a  volume  V  at  the  balance  temperature.    In  this  special 
case,  with  6n  ■      =  0,  eq.  (22)  becomes 

0  =  (Mj-pAVx)  -  (DT-pAV)  .  (Al) 
DT/V,  we  can  rewrite  (Al )  as 

mx"Pavx  ■  ^-paV 


=  DT(l-pA/pD)  +  paVx  =  M1  .  (A2) 

Now  we  ask  the  question,  "What  if  one  assumes  the  mass  of  D  is  equal  to  NQ,  its 
nominal  value,  and  the  density  of  D  is  equal  to  pR,  the  basis  density  of  the  appar- 
ent mass  scale  to  which  the  balance  weights  have  been  adjusted?"    This  assumption 
is  not  correct  because 

PR  '  PD 

and 

DT  f  HD. 


Since  pn  = 


or, 


Al 


Nevertheless,  the  incorrect  assumption  would  tell  us  that 

mJ  =  nd(i-pa/pr)  +  pavxem2  . 


(A3) 


We  can  next  ask  ourselves  "Under  what  conditions  does  M-|  =  M2?" 

A 

Recalling  that  Mn  is  made  equal  to  Nn,  the  nominal  dial  value,  we  subtract  eq. 
(A3)  from  (A2).  U  U 

MrM2  =  DT  (l-pA/pD)  -  Mg(l-pA/pR)    .  (A4) 

Also,  by  definition,  (see  eq.(  11)  in  the  main  text), 

=  DT  (1-P0/PD0)/(1-P0/PR0)  (A5) 

where  the  subscript  zeros  refer  to  the  standard  conditions  specified  in  the  defini- 
tion of  apparent  mass.    Substituting  (A5)  into  (A4), 

M2-M1  =  DT[(l-p0/pD0)(l-pA/pR)/(l-p0/pR0)  -  (l-pA/p0)]  • 

The  first  approximation  we  will  make  is  that  pRq  =  pR  and  pqq  =  pp.    That  is, 

the  temperature  is  sufficiently  close  to  20  °C  that  the  expansion  in  volume  of  the 
balance  weights  is  negligible.  Thus 

M2-M1  =  DT[(l-p0/pD)  (1-pa/pr)/(1-p0/pr)  -  0-pa/pd)]  • 
Next,  we  use  the  relation  developed  in  the  main  text: 


-1  2  3 

O-p0/pr)    =  l  +  Pq/Pr  +  (pq/pr)   +  (Po/pR^  + 


so  that 


M2-M1  =  DT[(l-p0/P[))(l-pA/pR)  (l+p0/pR+(p0/pR)2  +  (p0/pr)3  +  ...)  -  (1-pa/pd)] 

2  2 

Finally,  we  multiply  out  but  throw  away  all  terms  of  order  (pq/pr)  >  (pAp0^pR^' 

2  3  3 

Pq/pdpr,  and  smaller.    Since  pQ  ^   pa  <_  1.2  mg/cm   and  pD  ^  pR  >  7.5  g/cm  (the 

balance  manufacturer  has  seen  to  this),  then  neglecting  these  terms  leads  to  errors 

of  3  x  10"    percent  or  less. 

Thus,  we  now  have 

M2-M-|  =  DT  (p0/Pr  "  Pq/Pq  "  Pa/Pr  +  Pa/qD^ 
or,  more  simply 

M2-M1  =  DT  (p0-pA)  (l/pR-l/pD)  (A6) 

A2 


The  balance  manufacturer  has  selected  the  metal  for  his  weights  from  stock 
having  a  density  pd  sufficiently  close  to  pR  so  that,  near  sea  level,  M-.-M2  will  not 
exceed  the  balance  tolerance. 

In  recent  years,  balance  manufacturers  have  been  making  weights  of  stainless 
steel.    Equation  (A6)  shows  why  it  is,  therefore,  desirable  for  these  weights  to  be 
adjusted  on  the  8.0  basis. 

The  chief  virtue  of  the  apparent  mass  scheme  is  that,  even  though  different 
balance  manufacturers  may  make  their  weights  of  different  alloys,  a  user  can  be 
completely  unaware  of  these  subtleties  and  still  derive  reasonably  accurate  true 
mass  results  on  any  balance  so  long  as  the  correct  basis  density  is  used. 

Example: 

A  balance  has  20  g  of  built-in  dial  weights.    The  density  of  these  weights  is 
3 

actually  7.8  g/cm  .    The  weights  have  been  adjusted  so  that  the  dial  readings 
correspond  to  apparent  mass  on  the  brass  basis.    What  is  the  value  of  M2-M-|  at 

maximum  load  if  p^  =  1.16  mg/cm3?   What  if  the  laboratory  is  at  high  elevation  so 

that  pA  =  1.00  g/cm3? 

Answer: 

a)  M2-M1  =  20  g  (0.00120  -  0.00116)  [g^r  -  773]  =  -7  ug. 

b)  M2-M1  =  20  g  (0.00120  -  0.00100)  -  j^l  =  -37  ug. 

By  comparison,  the  calibration  uncertainty  for  high-quality  20  g-weights  as  measured 
by  NBS  is  about  10  yg;  the  tolerance  for  20-g,  Class  1,  metric  weights  is  74  ug 
[14]. 
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APPENDIX  B 


OUTLINE  OF  THE  DERIVATION  OF  p.,  THE  DENSITY  OF  AIR 

In  order  to  find  p^,  an  equation  of  state  involving  temperature,  humidity,  and 
barometric  pressure  has  been  developed.    We  will  follow  the  definitive  derivation  of 
F.  E.  Jones  [2]  and  briefly  describe  his  arguments  to  arrive  at  a  simplified  re- 
lation sufficiently  accurate  for  our  applications.    The  reader  is  urged  to  consult 
[2]  for  full  details. 

In  order  to  derive  the  density  for  a  mixture  of  two  gasses,  let  us  start  with 
the  ideal  gas  law. 


and 
where 


pV  =  nRT 


m 

n  =  M 


n  =  number  of  moles 

m  =  mass  of  gas 

M  =  molecular  weight  of  gas 

R  =  universal  gas  constant 

T  =  temperature  in  kelvins  =  (273.15  +  t)  for  t  in  degrees  celsius 


From  the  above,  we  can  rewrite 


and 


pV  =  jj|  RT  so  that  m  =  ^ 


Mm  PA  .  m 
—  =  —  =  —  since  pA  =  - 

RT     pV       p  V 


Specifically,  consider  air  consisting  of  dry  air  and  water  vapor.  For  the  dry  air, 
we  have 


mD  = 


PDVMD 


RT 

and  for  the  water  vapor,  we  have 

p  VM 
^w  w 
m  ,  =  — - —  . 

w  RT 

Since  the  total  density  is  given  by 

mD  +  mw 

pA  = 


Bl 


we  obtain 


PA  = 


PDVMD  Pwvmw 


/v 


RT 


RT 


RT 


Using  Dalton's  Law  for  partial  pressures 

P  =  P]  +  P2 


we  can  substitute  for 


so  that 


and  using 


PA  = 


1_ 
RT 
M 


PD  =  P  -  Pw 

(RMD  -  PwMD  + 
M, 


z  = 


M, 


M 

PA  =  IT  (P  +  PW(£  "  ]) 


RT 


At  this  point  it  is  important  to  consider  for  a  moment  the  correction  necessary  for 
a  gas  which  is  not  ideal.    To  do  so  we  rewrite  the  last  relation  as 


P  = 


PART 


1 


l+^(c-l) 

p 


For  an  ideal  gas  the  ratio 


pA 


RT 


=,=  Z 


1     +    ^  ( 

P 


-  1) 


B2 


has  to  be  equal  to  1 .  If  we  want  to  consider  non-ideal  gas  corrections,  then  we 
have  to  incorporate  Z  into  our  formalism.  Hence, 


P  = 


PARTZ 


1 


l+^(e-l) 

P 


where  Z  is  called  the  "compressibility  factor." 

Rewriting  this  relation  for  the  mixture  density,  we  now  have 

M 

-  [(P  +  (e  -  1)PW] 


RTZ 


Let  us  now  turn  to  the  dependence  of  this  relation  on  the  relative  humidity.  We 
define 

PW 

U  =  —  x  1 00 


where 

U    =  relative  humidity  in  percent 

p^l  =  effective  vapor  pressure  of  water  in  moist  air 

e£  =  effective  saturation  vapor  pressure  of  water  in  moist  air, 

Furthermore,  we  know  that 


where 


es  =  the  saturation  vapor  pressure  of  pure-phase  water, 


The  ratio  of  the  two  pressures  is  called  the  "enhancement  factor"  for  saturated 
water  vapor  and  is  given  by 


f  = 


so  that 


and  since 


es  =  es  *  f 


U  •  el 


100 
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we  have 


PW 


es  •  f 


100 


Substitution  into  the  mixture  density  relation  yields 


Pa  = 


RTZ 


[P  +  (e  -  1) 


U-e  -f 

— 
100 


where  U  must  be  given  in  %. 

This  relation  constitutes  the  final  formula.  The  parameters  P,  T,  and  U  must 
be  measured  by  the  user. 

Substituting  the  best  available  values  for  R  anf  Mw;  and  choosing  reasonable 
values  for  MD,  Z,  and  f,  which  are  approximated  as  constant  parameters,  we  have 

R  =  8.31441  J'ou1es    =  8,314.4  j°u1es 
mole-K  kmole-K 

MD  =  28.964  g/mole 

Z  =  0.9996 
f  =  1.0042 
Mw  =  18.0152  g/mole. 


Since 


£  = 


_  MW  _  18.0152 


MD  28.964 


we  obtain 


pA 


=  3  • 4848    (p  .  0.0037960  UeJ  x  10"3 
T  s 


mg_\  =  0.0034848    (p  _  QM3mQ  UeJ 


cm' 


where 


(t+273.15) 


t  is  in  °C 

P  is  in  pascals  (133.3224  Pa  =  1  mm  Hg)  . 


Converting  to  mm  Hg  pressure,  we  have 
pA 


ms\.    0.46460     (P  .  0.0037960  UeJ  . 
3      (t+273.15)  s 


cm 
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The  parameter  eg  has  been  determined  by  fitting  measured  data  between  288.15  K  and 
301.15  K.    The  relationship  developed  is 

-5315.56 

es  =  (1.3146  x  109)  eU+273'15)_mm  Hg. 

It  should  also  be  noted  that  the  value  for  the  enhancement  factor  used  above, i.e. 
f  =  1.0042  can  be  approximated  more  accurately  by 

f  =  1.00070  +  4.150  x  10"6  x  P  +  5.4  x  10"7t2  . 

In  the  above  final  equation,  the  relative  humidity  is  given  in  percent 

U  =  %  (e.g.  51.2,  43.7  etc.) 

and  the  pressure  P  in  mm  Hg. 

MD  (and,  therefore,  e)  depends  on  the  mixture  of  gases,  other  than  water  vapor, 

which  makes  up  ambient  air.    The  chief  variability  in  this  mixture  comes  from  CO2 

and  O2  levels.    These  are  assumed  perfectly  correlated—that  is,  CO2  levels  can  only 

increase  locally  at  the  expense  of  O2  and  vice  versa--as  in  processes  of  combustion, 
respiration  and  photosynthesis. 

Fortunately,  the  variation  has  little  effect  on  Mp.    In  eqs.  (20a)  and  (20b) 
(main  text)  we  have  assumed  the  ambient  level  of  CO2  typically  found  in  the  NBS  mass 
laboratories  (0.00042  mol  of  C02/mol  of  air).    A  100  percent  increase  in  this  level 
would  raise      by  less  than  0.02  percent. 

Z  is  slightly  dependent  on  barometric  pressure,  temperature  and  relative 
humidity  [2].    Between  19  °C  and  26  °C;  525  mm  Hg  and  825  mm  Hg;  and  0  and  100% 
R.H.,  Z  varies  by  less  than  0.03  percent. 
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APPENDIX  C 

PROPAGATION  OF  ERROR  THROUGH  THE  AIR  DENSITY  EQUATION 


Using  eq.  (20b)  (developed  in  Appendix  B),  we  have 

=    0- 46460  (p.o.o037960-U.ec) 
M     (t+273.15)  s 

For  ease  of  writing,  we  will  set 

P  =  PA 
a  =  0.46460 
b  =  273.15 
c  =  0.0037960 
e  =  es 

such  that 

p  =  — — (P  -  c-U-e)  . 
(t+b) 

We  want  to  determine  the  uncertainty  for  p  using 

dp  =  l£  dt  +  l£  dP  +  9£  du  . 
at         9P  9U 

at 


Using  (CI ) ,  we  find 


^=  a(P  -  c-U-e)  ^-  (— ) 
at  at  t+b 

t+b 


ap 


Using  (CI ) ,  we  find 


ap  _  a 
3P  t+b 


and  since 


p  ^ 


(t+b) 

CI 


where  we  ignored  the  second  term  in  (CI)  since  it  is  only  a  secondary  source  of 
error,  we  have 

.  (C4) 

8P  P 

3p_. 
9U 


From  (CI ) ,  we  obtain 

aP  acel) 

p  

t+b  t+b 

so  that 

9p  .  ace 
dU  (t+b) 

aP 

and  again  using  p 
we  have 

.  (C5) 

3U  P 

Next,  we  have  to  estimate  the  combined  uncertainties  due  to  t,  P,  and  U.    In  Section 
2.3.2  (main  text)  we  asserted  that  the  temperature  should  be  known  to  ±0.4  °C  so 
that 

dt  =  ±0.4  °C  . 

Similarly, 

dP  =  ±1.1  mm  Hg 

and 

dU  =  ±16%  . 

3 

All  these  values  assume  that  we  want  to  know  the  air  density  to  0.0017  mg/cm  . 
Substitution  of  all  terms  into  (C2)  yields 

dp  =  =2-  (±0.4)  +  £  (±1.1)  -  £e  p  (±16)  . 
t+b  P  P 

In  the  worst  case,  these  uncertainties  would  add  linearly.    However,  we  assume 
the  errors  are  uncorrected  so  that  a  better  estimate  of  the  total  uncertainty  is 
given  by   

-V(5*),*(5*)'-(S">' 


C2 


Jp2-(0.4)2  ~  P2(1.1)2  ~  C2e2p2(16)2 
*    (t+b)2  P2  P2 


16     +  TL21  +  C2e2256 


(t+b)2       P2  P2 


At  standard  conditions  we  have 


t  =  20  °C 

P  =  760  mm  Hg 


We  now  have 


"  (or 


x  L21    x  C2e2256 
<5p    =  Pi/  +   +   

(20+b)2     (-760)2  (760)2 


Substituting  p  =      and  replacing  b  and  c  with  their  numerical  values,  we  have 


^A=  J  0-16  +  1.21  +(0-0037960)2  ■  256-e2 
PA  (20+273. 15)2     (760)2  (760)2 


Also  es  at  20  °C  is  17.54  mm  Hg,  so  that  finally 


<$Pn 

—  -  0.0024 


pA 

or 

6pA 


0.24% 


C3 
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APPENDIX  D 

SOME  STATISTICAL  CONCEPTS 

We  start  out  with  an  infinitely  large  population  of  a  random  variable;  to  be 
practical,  let  us  say  a  very  large  number  of  similar  balance  observations.  We 
assume  that  this  population  is  normally  distributed  with  a  mean  of  u  and  a  variance 
of  c2. 

Let  us  take  a  random  sample  of  size  n  from  the  population.  Then  the  mean  of  the 
sample  is  defined  as 

X  =  -   E    X.  (Dl) 
n  i=l  1 

with  an  estimated  variance  for  a  single  observation  of 

n  ? 

(x.-xr 

s2  =  ^   (D2) 

n-1 

2  2 

As  n  grows  infinitely  large  %  -*■  u  and  s  ->  c  .  (This,  in  fact,  defines  what  we  mean 
by  u  and  a) . 

The  treatment  of  least  squares  data  is  more  complicated  than  analysis  of  repeat- 
ed measurements  of  the  same  quantity.    In  this  case,  average  values  of  several  dif- 
ferent masses  are  calculated  at  the  same  time  from  data  obtained  during  one  weighing 
design.    If  N  measurements  were  necessary  to  complete  the  design,  the  metrologist  has 
in  hand  the  results  of  these  N  mass  comparisons:    6-j ,  $2'  ...  6^.    The  least  squares 

analysis  provides  a  set  of  "best"  estimates:    6-j,  ^""^n*  as  we^  as  a  set  °f 

"best"  estimates:    mJ,  M^,...!^,  of  the  K  unknown  masses  calibrated  by  means  of  the 

design.    These  "best"  estimates  are  linear  combinations  of  the  measured  values 

6-j ,  fig,... wl'th  the  mass  of  the  known  standard  used  as  the  restraint.    The  linear 

combinations  are  uniquely  determined  by  the  choice  of  design  and  restraint.  Least 
squares  solutions  to  the  most  useful  weighing  designs  are  tabulated  in  ref.  [3]. 

The  standard  deviation  of  the  least  squares  fit  to  the  design  data  is  estimated 
from  the  formula^ 

N  2 

s2  =  ^   (D3) 

N-K+l 

where  the  denominator  is  given  by  the  degrees  of  freedom  (=  number  of  observations  - 
number  of  unknowns  +  number  of  restraints).    If  the  same  measurement  design  were 
repeated  m  times,  a  better  estimate  for  the  standard  deviation  of  the  process  could 
be  obtained: 


\east  squares  is  so  named  because  the  least  squares  solution:    6-j,  S^.-.S^, 


2 

minimizes  the  value  of  s    in  (D3). 


Dl 


s    =   (D4) 

m 

2  th  2  2 

where  s.  is  the  estimated  variance  of  the  i     run.    As  nh*»°,  s  -kj  .    This  defines 

2 

what  we  mean  by  o    in  the  case  of  least  squares.    Note:    the  number  of  degrees  of 

2 

freedom  (D.F.)  in  s    as  defined  by  (D4)  is  m(N-K+l).    Also,  this  a  is  the  same  o 
discussed  in  the  text,  the  "process  standard  deviation." 

Besides  a,  we  would  also  like  to  know  the  standard  deviation  of  each  individual 
mass,  M.,  computed  from  the  least  squares  analysis.    The  answer  involves  matrix 
algebraJand  depends  both  on  the  design  and  restraint  which  were  used,  as  well  as  on 
a.    All  of  the  commonly  used  cases  are^tabulated  in  ref.  [3].    In  general,  the 
standard  deviation  of  a  measured  mass  M.  is  some  value  £.0,  where  £.  is  a  number  that 

J  J  J 

depends  only  on  the  particular  design  and  restraint.    This  number  can  be  found  in 
ref.  [3];  0  is  defined  above. 

When  two  or  more  weights--j  and  p  for  example— are  used  in  combination,  their 

2       2  1/2 

combined  standard  deviation  is  not  simply  (£j  +  O      o-  but  is  often  somewhat 

greater.    This  is  because  data  taken  from  a  design  are  usually  correlated.    Both  the 
recipe  for  handling  weight  summations  and  tabulated  values  for  many  special  cases  can 
be  found  in  ref.  [3]. 

From  the  above  discussion,  it  should  be  evident  that  care  must  be  taken  in 
choosing  the  design  and  restraints  for  any  set  of  weights  to  be  calibrated.  Some 
choices  will  minimize  individual  I-  values  but  weight  summations  may  have  large 
uncertainties. 

2 

Assume  that  an  estimate  of  a   has  been  computed  from  m  designs  by  application  of 

2 

eq.  (D4).    Future  designs  will  produce  individual  values  s.  which  will  also  be 

2  J 
estimates  of  c   as  long  as  the  measurement  process  remains  stable.    An  F  statistic 

F  =  -A  (D5) 
s 

2  2 

can  be  computed  for  each  new  s..  Its  purpose  is  to  test  the  agreement  between  s. 

2  «J  2  " 

and  s    --or  more  precisely  to  test  whether  or  not  s-  comes  from  the  same  distribution 

2  J 
of  measurements  that  produced  s  . 

Given  a  large  number  of  F  statistics,  each  of  which  is  based  on  the  same  design, 

a  histogram  of  these  F  statistics  will  fall  very  nearly  on  the  universal  curve  known 

as  the  F  distribution.    The  theoretical  calculation  of  this  curve  depends  only  on  the 

2  2 
degrees  of  freedom  in  each  s.,  i.e.  N-K-l,  and  on  the  degrees  of  freedom  in  s  ,  i.e. 

J 

m(N-K-l).  The  distribution  is  scaled  so  that  the  area  under  the  curve  is  equal  to  1, 
making  it  possible,  for  example,  to  find  the  point      such  that  99  percent  of  all  F 

values  will  be  less  than  F^.    As  a  consequence,  the  next  time  that  we  carry  out  the 

weighing  design,  we  expect  that  there  is  a  99  percent  chance  that  the  F  statistic 
computed  from  that  particular  design  will  be  less  than  the  percent  point  F^.    If,  in 

fact,  it  turns  out  that  the  F  statistic  is  greater  than  the  percent  point  F^,  the 

D2 


precision  for  that  design  is  poorer  than  is  expected,  and  the  design  should  be 
repeated. 

The  percent  point  Ft  is  often  represented  in  tables  as  F(v-i ,  v^J-a)  because  it 
depends  on: 


v-j  =  the  degrees  of  freedom  in  s.  which  is  N-K+l 


V2  =  the  degrees  of  freedom  in  s^  which  is  m(N-K+l) 

a   =  the  significance  level  such  as  a  =  0.01  (i.e.  99  percent) 

For  mass  calibrations  at  NBS,  the  significance  level  is  chosen  to  be  a  =  0.01, 
and  because  the  measurement  process  has  been  tracked  for  a  long  time  resulting  in  an 
2 

estimate  of  a  that  has  a  ^ery  large  number  of  degrees  of  freedom,  the  value  that  is 
used  for  the  test  is 

Ft  =  F(vrcos  0.99) 
In  this  special  case,  Fx  can  be  wel 1 -approximated  by*: 


Ft  =  f  1   -  +  2.32635 


for 


and 


for 


9(N-K+1)  \  9(N-K+1) 


v]  >  2 


Ft  =  6.64 


=  1 


Tables  and  detailed  discussion  on  the  F  value  and  t  ratio  can  be  found  in  ref.  [13]. 


* 

A  derivation  of  this  result  is  well  beyond  the  scope  of  this  work.    The  interested 
reader  may  wish  to  consult:    Paulson,  Edward.    An  approximate  normalization  of  the 
analysis  of  variance  distribution.    Annals  of  Mathematical  Statistics.    13:  223-235; 
1942. 
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